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ABSTRACT 


The  influence  of  strakes  placed  along  vortical  flows 
has  been  determined  in  order  to  define  ways  of  drag  reduction 
on  wake-like  flow  fields  present  on  aircraft  afterbodies. 

Strakes  of  different  aspect  ratio  and  thickness  were 
tested  at  different  angles  of  incidence  with  respect  to  a  vortex 
path  in  order  to  evaluate  the  best  configuration  which  can 
reduce,  or  minimise,  the  cross  flow  kinetic  energy,  thus  the 
vortex  drag  of  a  vortical  flow. 

The  flow  field  in  the  vicinity  of  the  boattail  of  a 
scale  cargo  aircraft  was  studied.  The  aircraft  was  tested 
with  and  without  strakes  placed  on  its  afterbody,  the  results 
confirm  that  when  strakes  are  properly  placed  on  the  boattail, 
a  reduction  in  drag  is  achievable. 
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Wind  tunnel  and  flight  test  results  have  shown  that  when 
strakes  are  mounted  in  the  afterbody  part  of  an  aircraft  (or  so 
called  "boattail")  a  reduction  in  the  overall  drag  may  be 
obtained. 

Previous  researchers  have  investigated  the  flow  on 
afterbodies  of  representati ve  fuselages  of  high  upswept  angle. 
Different  devices  such  as  strakes  and  flow  deflectors  have  been 
fitted  to  such  configurations  and  qualitative  conclusions  from 
visualization  experiments  were  drawn.  In  some  cases  the  effect 
of  applying  suction  along  the  boattail  area  was  investigated 
and  the  effects  on  the  afterbody  Induced  wake  was  as  well  deter¬ 
mined  qualitatively. 

Although  helpful,  no  previous  work  done  on  aircraft 
afterbodies  could  define  exactly  what  is  the  effect,  and  its 
reason,  when  a  strake  is  placed  along  a  wcke-like  flow,  or 
vortical  flow,  like  the  ones  Induced  by  aircraft  with  a  large 
unswept  boattail  angle. 

This  prompted  the  desire  to  carry  out  a  series  of 
experiments  at  the  von  Karman  Institute  in  order  to  define  the 
mechanisms  of  action  of  a  strake  in  the  presence  of  vortical 
flow.s  and  it  application  on  aircrafts  as  a  way  of  drag  reduction. 
This  repor  -ibes  such  results. 

1  ’  . _ .earch  steps  carried  out 

In  order  to  carry  out  our  proposed  projects,  several 
steps  were  followed. 


First,  we  investigated  the  flow  field  existent  in  an 
aircraft  where  the  fusel age-boattai 1  combination  induces  a  vortex 


Secondly,  using  a  delta  wing  placed  at  a  high  angle  of 
attack,  the  flow  field  downstream  of  its  trailing  edge  was 
studied.  This  was  done  since  previous  research  (Ref.  16)  has 
shown  that  the  vortex  field  induced  by  aircraft  afterbodies  is 
similar  to  the  downstream  vortex  induced  by  separation  on  lift¬ 
ing  slender  delta  wings  at  high  angles  of  attack. 

By  placing  strakes  of  different  configurations  along 
the  vortex  path  downstream  of  the  delta  wing  trailing  edge,  and 
scanning  the  flow  behind  the  strakes  we  determined  the  effect 
of  strakes  in  a  vortical  flow,  and  most  important,  defined  ways 
of  drag  reduction  by  defining  which  strake  configuration  mini¬ 
mizes  the  vortex  cross  flow  kinetic  energy. 

The  above  information  served  us  to  know  how  strakes 
should  be  placed  on  aircraft  afterbodies.  This  was  done  and 
the  results  were  validated  by  measuring  the  static  pressures 
along  the  aircraft  afterbody  surface. 


2.  THEORY 


2 . 1  Determination  of  vortex  drag 

The  analysis  proposed  by  Maskell  (Ref.  5)  which  is  an 
extension  of  the  theory  due  to  Betz  to  calculate  the  drag  of  a 
wing  of  finite  span  from  a  wake  survey  in  one  plane  downstream 
of  the  model  was  adopted,  in  order  to  define  what  parameters  on 
the  wake-drag  equation  could  be  minimized  so  that  its  drag 
could  be  reduced. 


Consider  the  flow  generated  by  a  lifting  wing  of 
finitespan  in  a  closed  wind  tunnel,  where  u,v,w  are  the  x,y,z 
components  of  velocity,  P$  the  static  pressure,  Hy  the  total 
head  and  D  the  total  drag  of  the  wing.  If  the  flow  is  assumed 
steady,  incompressible  and  irrotational  outside  boundary  layers 
and  wakes  the  conservation  of  momentum  equation  within  a  control 
surface  bounded  upstream  and  downstream  of  the  wake  requires 
that 
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(1) 


This  equation  does  not  take  into  account  the  effect  of  pressure 
gradients  between  the  control  surfaces. 


The  total  head  formula  is,  by  definition 


H  =  P  +  -  p  (u2+v2+w2)  (inside  wake) 


and  outside  the  boundary  layer  and  wake  of  the  wing 


Ht  =  P  +  i  pU2 
T~  00  2 


(2) 


(3) 


where  Uw,  Pu  are  the  velocity  and  static  pressure  of  the  undis¬ 
turbed  stream. 


Substituting  equations  (2)  and  (3)  into  (1)  for  the 
two  control  planes  respectively  we  obtain  : 

f  f  1  C  2  2 

D]  =  jj  (Ht  -HTz)  dS  +  -  J  J  ( u ! _  u2  )  dS 
M 

♦  1  P  ft  f(V2+W2)-(v?+W?)]  dS  (4) 

2  J  J  l  J 

f 

where  denotes  the  integral  over  the  wake  region. 

.  . 

W 

Equation  4  implies  that  the  drag  of  a  wake  is  propor¬ 
tional  to  the  sum  of  the  cross  flow  velocity  components  squared 
at  the  control  surface  downstream  of  the  wake.  Also,  as  the 
value  of  the  U2  velocity  component  approaches  the  undisturbed 
stream  velocity,  the  second  term  of  the  RHS  equation  (4)  is 
kept  to  a  minimum. 

• 

Note  that  Dx  does  not  contain  all  the  terms  found  in 
a  detailed  vortex  drag  equation.  Reference  14  details  it 
including  wall  effects  and  blockage  terms. 


All  measurements  were  acquired  by  means  of  the  data 
acquisition  system  which  can  be  seen  in  figure  la.  It  consists 
of  a  floppy  disk  driver  linked  to  a  Commodore  PET  computer 
which  acquires  a  signal  at  a  predetermined  lapse  by  means  of  a 
stepping  motor-scani val ve  combination.  The  data  is  digitized, 
stored  in  a  digital  cassette  and  reduced  by  the  V K I  central 
computing  system.  The  data  is  reduced  as  well  by  the  PET  com¬ 
puter  and  some  of  the  results  are  prompted  during  testing  in 
order  to  check  the  validity  of  the  data  acquired.  All  pressures 
were  checked  against  a  Betz  manometer.  Figure  lc  shows  the  data 
acquisition  schematic  used  for  our  experiments.  Depending  on 
whether  it  was  surface  static  pressures  or  velocity  field  plots 
that  we  searched  for,  pressure  taps  or  a  five  hole  probe  were 
connected  into  a  multiple  portion  scan ivalve.  A  Validyne 
MP45-1  variable  reluctance  transducer  measured  the  pressure 
difference  with  respect  to  the  tunnel  static  pressure.  The 
voltage  signal  reached  an  amplifier  of  sufficient  gain  to  opti¬ 
mize  the  accuracy  of  measurements  :  an  analog  to  digital  conver¬ 
ter  was  placed  before  the  PET  Commodore  computer. 

Figures  2a  and  2b  show  the  five  hole  probe  used  for 
our  measurements.  It  was  supported  by  a  traversing  mechanism 
fixed  along  one  of  the  straight  walls  of  the  test  section  and 
flow  field  surveys  were  made  by  displacing  the  probe  along 
surface  planes  perpendicular  to  the  test  section.  The  stem 
which  encloses  the  probes  was  aligned  with  the  longitudinal 
axes  of  the  wind  tunnel.  A  detailed  view  of  the  five  hole 
probe  can  be  seen  in  figure  3  and  the  support  system  is  shown 
in  figure  4. 

The  calibration  curve  for  the  pressure  transducer  is 
shown  in  figure  5  showing  the  linearity  of  the  instrument;  its 
deviation  is  0.5%. 


Appendix  I  includes  the  data  acquisition  program 
implemented  in  the  PET  computer  which  was  used  for  our  tests. 

3.2  Measurements  using  a  five  hole  probe 

Flow  field  measurements  were  done  by  means  of  a  five 
hole  probe  (Figs.  2a, 3)  in  order  to  obtain  a  quantitative  des¬ 
cription  of  the  three  dimensional  flow  involved  in  our  case. 

This  is  of  prime  importance  since  we  deal  with  vortex  flows 
where  limiting  streamline  pattern,  smoke  injection  and  pressure 
plots  at  the  surface  are  not  sufficient  to  understand  the  flow 
behaviour . 

The  five  hole  probe  method  of  measurement  was  used 
since  it  can  be  set  in  a  non-nulling  mode  of  operation  (i.e., 
the  probe  is  inserted  into  the  flow  in  a  fixed  attitude),  and 
physical  quantities  describing  the  flow  can  be  deduced  from 
the  pressures  recorded  by  the  probe. 

Small  measurement  errors  in  vortical  flows  are  expected 
due  to  the  finite  probe  size  in  flow  gradients,  unsteady  effects 
and  probe  response  time. 

Disturbances  produced  by  the  probe  own  geometry  are 
accounted  for  by  calibrating  it  in  a  free  stream  of  known  speed 
and  direction.  The  calibration  speed  was  the  same  as  test  con¬ 
ditions  with  pitch  angle  limits  higher  than  the  expected  in 
vortex  flows. 


3.2.1  £a l^i brati_on  of_the_fi^ve  hoj[e_probe 

This  type  of  probe  enables  the  total  local  pressure 
and  velocity  vector  components  to  be  deduced  from  the  pressures 
sensed  by  the  five  tubes  open  to  the  flow  in  a  configuration 
shown  in  figure  3. 


The  calibration  was  done  following  a  combination  of 
the  methods  described  by  Wickens  (Ref.  10)  and  Meineke  (Ref.  7) 
The  parameters  listed  below  are  deduced  from  the  pressures 
measured  by  the  probe.  These  are  : 


PD 


P5-P 


dynamic  pressure  parameter 


,  V  2 


CP  = 


PT-P5 


-  PV2 


total  pressure  parameter 


PHI 


P2-P4 


PI  - P3 


roll  parameter 


BIR 


n/(P2-P4)  +  ( PI  -P3  )2 
P5-P 


misalignment  parameter 


where  P1-P4  are  the  pressures  at  the  probe  defined  in  figure  6B 
These  relations  are  sufficient  to  obtain  a  calibrated  probe  if 
the  flow  direction  is  given  by  the  pitch  angle  e,  and  the  roll 
angle  $  normal  to  the  probe  head.  The  actual  roll  flow  angle 

-i  P2-P4 


can  then  be  deduced  from  the  relation  tan" 


P1-P3 


The  pitch  and  roll  angles  and  the  flow  direction 
quadrants  are  defined  according  to  the  convention  shown  in 
figure  6b.  The  calibration  duct  used  for  our  probe  is  shown  in 
figure  6a.  By  keeping  its  velocity  constant  and  varying  the 
pitch  and  roll  angles  covered  by  the  different  quadrants,  the 
calibration  data  shown  in  figures  7-11  was  finally  obtained. 
These  curves  were  approximated  by  a  N  degree  polynomial  (its 
degree  was  chosen  depending  on  the  data  scatter)  and  input 
in  the  data  reduction  program.  Each  of  the  plots  in  figures 
7-11  show  the  experimental  data  curves  and  points  (dashed  line) 
and  the  polynomial  curve  (solid  line)  which  was  used  to  approxi 
mate  the  data.  The  polynomial  formulas  used  can  be  read  from 
Table  1. 
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Figure  10  shows  the  agreement  between  the  calculated 
roll  angle  and  the  one  set  in  the  calibration!  duct 

The  dynamic  pressure  parameter  variation  with 
respect  to  the  roll  angle  can  be  seen  in  figure  11.  Although 
not  as  symmetric  as  Wickers'  results,  its  cyclic  variation 
with  quadrants  variations  is  well  defined.  / 

3.2.2  Calculation  of_total_pressure  and 
f 1 ow_vel ocity  components 

During  tests,  the  central  and  circumferential  pres¬ 
sures  measured  from  the  five  hole  probe  were  recorded.  These 
results  were  used  to  calculate  4> ,  deduce  the  quadrant  from 
which  the  flow  is  coming,  and  the  value  for  BIR.  Knowing  BIR, 
| 0 |  could  be  found  from  the  calibration  curves. 

Secondly,  the  parameter  PD  could  be  obtained  from 
the  calibration  polynomials.  This  allowed  us  to  deduce  the 
value  of  the  local  velocity  vector  using  the  following 
relation  : 


V  = 


The  parameter  CP  was  then  found  from  the  calibration 
plots.  Finally,  the  total  local  pressure  could  be  obtained  by 

using  PT  =  P5  +  -  pV2CP. 

2 

The  probe  axis  system  and  the  Cartesian  system  are 
shown  in  figure  6b.  The  transformation  formulas  used  were  : 


s  i  n  |  ct  I  =  s  l  n  |  e  1  x  c  o  s  j  4>  | 


cos  |  y  | 


_ cos | e _ 

•v/l  -si  n2  |  e  j  xcos2  |  $  | 


and  vice  versa 


tan |  <j>  |  =  —  n-Ly~- 
tan|o| 

cos | 9 |  =  cos |a |  •  cos | Y | 

Velocity  component  signs  were  given  by  knowing  the 
flow  quadrant.  The  quadrant  numbering  convention  was  defined 


as  follows  : 

Quadrant  I:  a  >  0°  Y  >  0° 

PI  >  P3 ,  P4  >  P2 

Quadrant  II:  a  >  0°  Y  <  0° 

PI  >  P3 ,  P2  >  P4 

Quadrant  III:  a  <  0°  Y  <  0° 

P3  >  PI,  P2  >  P4 

Quadrant  IV:  a  <  0°  Y  <  0° 

P3  >  PI,  P4  >  P2 


Finally,  the  flow  velocity  components  were  calculated 
using  the  following  relations  : 

u  =  V  x  cosa  x  cosy 

v  =  V  x  cosa  X  sinY 


w  =  V  x  sina 


The  L-2A  is  an  open  return  closed  test  section  wind 
tunnel  as  shown  in  figure  12.  The  speeds  at  which  all  our 
measurements  were  taken  was  at  29  m/s  (50  mm  H20)).  In  order 
to  ensure  that  the  cross  flow  velocity  irregularities  present 
in  the  clean  tunnel  would  not  influence  the  accuracy  of  our 
experiments,  two  flow  surveys  were  made  at  different  test  sec¬ 
tion  stations;  measurements  of  turbulence  were  carried  out  as 
well.  Figure  13  shows  the  cross  sectional  view  of  the  L-2A. 

The  solid  line  bounds  the  area  where  the  flow  surveys  were 
made  and  the  dashed  line  is  the  boundary  of  the  turbulence  level 
measurements.  All  flow  field  and  turbulent  maps  are  shown 
looking  upstream  of  the  test  section. 

3.3.1  £low_qy.al_i  ty 

Figure  14  shows  a  flow  field  map  of  one  of  the  test 
section  stations  surveyed.  The  magnitude  of  the  cross  flow 
velocities  is  of  the  order  of  0.50  m/s,  which  gives  an  average 
flow  angularity  on  the  tunnel  of  about  1.0  degree.  The  rota¬ 
tional  motion  seen  in  figure  14  is  thought  to  be  inherent  of 
tunnels  of  such  classes.  The  magnitude  of  the  cross  flow  velo¬ 
cities  is  small  compared  to  the  results  of  our  experiments, 
therefore  no  large  margin  of  error  was  expected.  Isolines  of 
total  pressure  are  also  shown  across  the  test  section.  In 
case  one  is  to  compare  the  magnitude  of  the  cross  flow  vectors 
by  looking  at  the  length  of  the  arrow-vector  itself,  the  reader 
should  note  that  the  cross  flow  vector  plot  of  figure  14  is  not 
plotted  in  the  same  proportion  as  the  rest  of  the  cross  flow 
vector  plots  existing  in  this  report. 

3.3.2  Tu  rbu  1  encejnea  s_ur  eme  n  t£ 

Figures  15a  and  15b  show  turbulence  levels  and  iso¬ 
lines  found  at  the  same  station  as  the  cross  flow  survey  was 
done.  The  average  percent  of  turbulence  found  at  the  center 
of  the  test  section  was  of  about  0.50%. 


Figure  16  shows  isolines  of  P  ,  found  in  the  L-2A 
test  section. 

Figures  17a  and  17b  show  a  schematic  from  the  models 
tested  during  our  research  as  they  were  positioned  in  the 
tunnel  with  the  boundaries  of  flow  quality  survey  and  turbulence. 

Geometric  blockage  calculations  indicate  that  for  the 
models  tested  their  value  was  of  3%. 

3 . 4  The  aircraft  model 

In  order  to  study  the  flow  field  in  the  vicinity  of 
an  aircraft  afterbody,  a  1/72  scale  model  of  a  Lockheed  Hercules 
C-130  E  was  used.  Figures  18a  and  18b  show  the  aircraft  tested, 
the  upswept  angle  at  the  fusel age-boattai 1  junction  is  of  approxi 
mately  25°  and  reaches  values  as  high  as  30°  at  the  end  of  the 
boattail.  Forty  eight  pressure  taps  were  installed  in  a  circum¬ 
ferential  manner  along  the  fuselage  boattail  on  four  different 
fuselage  stations  (12  orifices  along  each  station)  and  their 
values  measured  by  means  of  tubing  installed  through  the  air¬ 
craft  wing  (Fig.  57).  All  tests  on  the  aircraft  were  done  at 
zero  degrees  fuselage  angle  of  attack.  The  Reynolds  number 
of  tests  based  on  fuselage  diameter  was  of  1.8  x  105.  Transi¬ 
tion  strips  which  consisted  of  sand  of  height  dimensions  vary¬ 
ing  between  .420  and  .590  mm  were  placed  along  the  beginning 
of  the  fuselage  and  on  the  upper  surface  of  the  wing.  The 
first  one  was  placed  in  order  to  ensure  having  turbulent  flow 
at  the  afterbody;  the  transition  strip  over  the  wing  was  placed 
to  create  the  correct  downwash  angle  of  the  wing  along  the 
boattai 1  . 

Figure  19  shows  the  experimental  set  up  as  the  model 
was  installed  in  the  test  section.  For  measuring  the  static 
pressures  on  the  boattail,  the  48  port  scanivalve  placed  on 
the  lower  left  corner  was  used.  The  traversing  mechanism 
seen  on  the  right  was  used  to  displace  the  five  hole  probe. 
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A  three  view  picture  of  a  C-130  model  is  seen  in 
figure  20.  The  side  view  denotes  the  different  fuselage 
stations  where  the  flow  field  was  surveyed.  Fuselage  station 
No.  1  is  located  at  the  fuselage-afterbody  junction  and  No.  5 
is  at  the  end  of  the  boattail  whereas  No.  6  is  located  one 
fuselage  diameter  downstream  of  No.  5. 

Note  that  fuselage  station  3  of  the  aircraft  corres¬ 
ponds  to  the  same  wind  tunnel  longitudinal  station  at  which 
the  flow  quality  survey  and  turbulence  measurements  were  made. 

Pressure  taps  were  installed  along  fuselage  stations 
Nos.  1,2, 3, 4  and  their  positioning  can  be  seen  in  figure  57. 

3.4.1  0i'2  smoke  v.i  sual  i_zat£on 

By  inserting  oil  smoke  upstream  of  the  aircraft  the 
boattail  induced  wake  could  be  visualized.  Such  is  the  cace  of 
figure  21  which  was  taken  at  fuselage  station  No.  3.  Even 
though  photographic  quality  is  not  the  best  since  the  pictures 
were  taken  with  a  camera  outside  the  test  section  focused  on  a 
mirror  inside,  the  wake  could  clearly  be  seen  which  also  gave 
us  an  estimation  of  the  measurement  domain  of  our  flow  field 
surveys.  The  smoke  flow  visualization  results  agree  qualita¬ 
tively  with  the  flow  field  surveys  using  the  five  hole  probe. 

3.4.2  £1  ow_f  i_e  l_d_s  u  r  v  ey  s 

In  the  case  of  the  C-130  where  there  is  a  high  surface 
curvature  at  the  junction  of  the  fuselage  side  panels  and  the 
afterbody  lower  area,  the  three  dimensional  boundary  layer  flow 
seems  to  leave  the  surface  near  this  junction  and  roll  up  into 
trailing  vortices  which  travel  along  the  boattail  as  they  are 
swept  upwards.  Such  is  the  indication  of  the  flow  field  surveys 
starting  at  fuselage  stations  No.  1  to  No.  5.  A  complete  wake 
can  be  deduced  to  exist  at  the  end  of  the  boattail  area  (No.  5) 


and  this  is  validated  when  the  flow  field  is  surveyed  downstream 
of  the  aircraft  (F.S.  No.  6).  The  measurement  domain  could  not 
go  higher  since  the  horizontal  tail  interfered  with  the  probe 
stem . 

Isolines  of  Pr  indicate  two  areas  of  pressure  loss 
at  the  end  of  the  boattail  as  indicated  by  figure  26b.  The 
first  is  along  the  side  of  the  boattail  and  the  second  of 
slightly  less  loss,  is  a  bluff  wake  which  leaves  the  fuselage 
right  behind  it.  When  isolines  of  cross  flow  kinetic  energy 
are  plotted  non  di mens i ona 1 i zed  with  respect  to  the  free  stream 
velocity  squared,  the  location  of  the  wakes  is  more  easily 
deduced,  as  shown  in  figure  26c. 

The  flow  field  survey  downstream  of  the  boattail 
(Fig.  27a)  shows  the  final  geometry  of  the  wake  as  it  leaves 
the  aircraft  in  conjunction  with  other  small  wakes  like  the 
one  near  the  vertical  tail.  Isolines  of  Pr  (Fig.  27b)  denote 
two  areas  of  large  pressure  loss  which  probably  are  the  boat- 
tail  induced  trailing  vortex,  and  the  bluff  wake  which  leaves 
the  aircraft  through  the  corner  of  the  end  of  the  boattail. 

Figure  27c  agrees  with  the  previous  figure  by  locat¬ 
ing  the  position  of  the  boattail  induced  vortex  with  isolines 
of  cross  flow  kinetic  energy. 

3.4.3  L i mi  ti  ng_s t^reaml  i_nes_vi_sua l_i  z.a t_i  on 

Limiting  streamlines  visualizations  was  done  by  mixing 
talc,  oleic  acid,  titanium  dioxide,  gasoline  and  oil  and  placing 
small  dots  on  the  afterbody  surface.  Figure  28a  shows  the 
top  view  of  the  aircraft,  where  limiting  streamlines  travel¬ 
ling  close  together  tend  to  separate  in  different  directions. 
Figure  28b  gives  an  idea  of  the  embracing  of  the  flow  over  the 
boattail  as  well  as  the  downwash  angle  induced  by  the  wing. 

Two  separation  lines  can  be  seen,  one  at  the  s i de-panel -bottom 
junction  and  the  other  which  is  also  seen  in  figure  28c, 


correspondi ng  to  the  bluff  wake  that  leaves  the  boattail 
through  the  afterbody  corner.  Note  on  figure  28c  the  result 
of  stagnated  flow  present  on  the  upper  end  of  the  boattail  ; 
some  of  the  mixture  used  to  visualize  the  limiting  streamlines 
did  not  move. 

3 . 5  The  delta  wing 

A  delta  wing  of  aspect  ratio  =  1.0  (Fig.  29)  placed 
at  20°  angle  of  attack  was  positioned  on  a  flat  plate  and  tested 
i n  the  wi nd  tunnel . 

3.5.1  F 1  o w  f ^ej_d_su rv^ey s_downs t r earn  of  it 

Surveys  were  made  in  order  to  determine  the  position 
of  the  vortex  path  as  it  travels  downstream.  Figures  30a, b,c 
show  the  flow  field,  Pp  and  CFKE  isolines  of  a  plane  £  =  0.625 
downstream  of  the  trailing  edge.  Figures  31a, b  show  the  cross 
flow  field  and  CFKE  values  of  a  plane  £  =  1.875  downstream  of 
the  trailing  edge. 

Results  of  the  flow  field  survey  at  5  =  2.5  are  shown 
in  figures  32a, b  and  c.  The  control  plane  has  been  increased 
compared  to  the  ones  shown  in  figures  30  and  31. 

The  i  =  2.5  station  is  the  basis  for  our  analysis 
when  we  will  compare  these  results  against  the  same  control 
plane  but  having  strakes  placed  before  it. 

3.5.2  Vortex  core  travel  path 


Figure  33  shows  the  travel  path  that  the  delta  wing 
induced  vortex  follows  as  it  leaves  the  trailing  edge.  Note 
that  the  helix  makes  an  angle  with  respect  to  the  momentum 
flow  direction  and  that  as  it  travels  downstream  it  tends  to 
curve  referenced  to  the  tunnel  axis.  The  height  of  the  vortex 
core  is  approximately  constant  as  it  moves  downstream. 


3.6.1  Thin  strakes 


Different  strake  configurations  were  placed  downstream 
of  the  delta  wing,  and  pivoted  at  £  =  1.25,  in  order  to  change 
the  angle  of  incidence  with  respect  to  the  incoming  vortex. 

Figure  34  shows  how  strakes  were  placed  on  the  flat 
plate  and  their  position  with  respect  to  the  wing.  The  flow 
field  downstream  of  the  strake  was  studied  (at  the  £  =  2.50 
plane) . 


Measurements  of  static  pressures  on  the  strake 
surfaces  were  also  done. 

Figure  35  shows  the  different  strake  configurations 
that  were  tested.  Water  tunnel  flow  visualization  studies  were 
made  by  placing  a  delta  wing  of  AR  = 1  at  20°  angle  of  attack 
and  observing  the  vortex  path  and  the  field  downstream  of  the 
strakes  qualitatively.  Figure  36  shows  the  set  up  used.  The 
vortex  formed  by  the  delta  wing  leading  edge  separation  as  it 
travels  downstream  was  observed  and  strakes  were  positioned 
along  the  path. 

The  height  of  the  strake  was  chosen  so  that  it  would 

reach  three  different  levels  of  H  (as  defined  in  figure  37). 

If  we  look  upstream  on  the  Pr  isolines  plot  of  figure  32b,  wake 

dimension  D  is  the  distance  from  the  vortex  core  to  the  P„.  =  1.1 

r 

isoline  (the  largest).  Therefore  H  =  0  when  the  strake  reaches 
P-r  =  1.1  on  the  lower  side  and  1.0  when  it  reaches  the  highest 
point  along  the  curve,  H  =  0,5  when  the  height  of  the  strake 
reaches  the  vortex  core.  The  angle  of  attack  sign  can  be  de¬ 
fined  in  the  conventional  way.  It  is  positive  when  the  leading 
edge  of  the  strake  is  moved  right  and  negative  when  it  moves 
left.  The  zero  alpha  position  was  when  the  strake  is  aligned 
and  in  the  center  of  the  vortex  path. 
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Figure  38  shows  that  when  the  strake  H  =  0  is  placed 
at  a  negative  angle,  as  the  helix  moves  along  the  strake  the 
lower  part  of  the  vortex  which  cannot  jump  to  the  other  side 
continues  having  a  vortical  type  motion,  but  the  flow  coming 
out  of  the  forward  facing  side  is  streamlined  as  a  result  of 
not  having  rotational  energy  or  vorticity  contributed  by  the 
flow  on  the  opposite  side  of  the  strake. 

The  same  sort  of  effect  is  present  when  ar  =  0, 
although  the  streamlined  flow  has  changed  positions  and  the 
main  vortex  core  has  been  displaced. 

For  positive  incidence,  the  appearance  of  an  opposite 
direction  vortex  field  from  the  main  appears  in  the  lower  right 
corner  with  small  values  for  its  cross  velocities.  For  negative 
incidences,  the  values  of  the  velocities  are  higher  than  for 
angles  of  zero  or  above. 

The  same  trend  is  followed  by  strakes  of  ratio  H  =  0.5 
(Figs.  41-43)  although  the  magnitude  of  their  cross  flow  veloci¬ 
ties  is  smaller  than  for  the  H  =  0  case,  especially  when  ar  =  0, 
which  produces  a  streamlining  effect  clearly  noticed. 

At  positive  angles  of  attack,  a  double  vortex  is 
definitely  present  which  reduces  the  magnitude  of  opposite 
direction  velocities  of  the  velocity  vector  integral,  thus 
reducing  the  vortex  energy  loss. 

Water  tunnel  visualization  for  H  =  0.5  strake  is 
shown  in  figure  44  which  confirms  the  result  shown  by  the  cross 
flow  velocity  components  plots  :  an  existing  vortex  is  still 
present  at  =  -15  and  is  streamlined  as  ar  =  0;  at  positive 
angles  of  attack,  part  of  the  vortex  flow  is  redirected  due 
to  the  influence  of  the  strake  presence  and  still  maintains 
vortex  motion,  and  the  other  plot  "jumps  over"  (Fig.  44b)  the 
strake  with  its  core  finding  a  new  equilibrium  position  near 
the  lower  corner  of  the  control  plane. 


V  v. 


Similar  observations  are  found  for  strakes  of  H  =  1.0 
(Figs.  45-47)  with  the  difference  that  the  vel oci ty  components 
are  higher  than  H  =  0.5  and  H  =  0. expect  for  ar  =  -15°.  Strakes 
of  H  =  1.0  have  a  wider  area  of  pressure  loss  where  a  vortex 
core  is  present.  This  can  be  observed  on  the  Pr  isoline  plots. 

Water  tunnel  visualization  for  strakes  of  different 
H  and  constant  a  is  seen  in  figure  48.  For  the  case  of  strake 
placed  in  the  lower  edge  of  the  vortex,  the  wing  induced  vortex 
does  not  get  influenced  as  much  as  the  other  two  cases  do. 

For  the  H  =  1.0  case  the  vortex  impinges  on  the  forward  facing 
surface  and  does  create  a  more  disorganized  flow  than  the 
H  =  0.5  case. 

3.6.2  Th i^c k_  strake s^ 

In  order  tc  evaluate  the  effect  of  thickness  on  the 
cross  flow  kinetic  energy,  strakes  of  H  =  0.5  ratio  were  tested 
at  different  angles  of  attack.  The  resultant  flow  field 
agrees  in  trend  and  in  approximate  magnitude  with  the  thin 
strakes  placed  at  similar  incidences.  Results  show  that  at 
ar  =  0  the  flow  field  is  much  more  streamlined  than  ap  =  -15 
and  that  for  «r  =  15®  a  double  vortex  formation  of  opposite 
velocity  direction  occurs. 

3.7  Pressure  measurements  on  strakes 


Figure  52  shows  the  three  different  thick  strake 
configurations  that  were  tested  and  where  static  pressure 
measurements  were  made  on  each  surface.  The  reason  for  these 
measurements  was  to  find  out  which  strake  configuration  bends 
providing  a  reduction  in  CFKE ,  would  at  the  same  time  create 
lift  or  favourable  force  along  its  axis. 

Figures  53,54  and  55  show  pressure  distribution  for 
H  =  0,  0.5,  1.0  strakes  at  three  different  angles  of  attack. 
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The  results  are  similar  for  the  three  strake  configu¬ 
rations.  This  implies  that  the  height  of  the  strakes  is 
irrelevant  for  its  pressure  sensed.  When  strakes  are  placed 
at  angles  of  attack  with  respect  to  the  incoming  vortex  of  0° 
or  less,  the  net  result  is  an  increase  in  lift,  which  as  a 
result  provides  an  increase  in  thrust  along  the  strake  axis. 

Calculations  of  lift  integrals  are  not  possible  since 
the  amount  of  pressure  drop  was  not  considered  enuugh.  The 
results,  though,  give  an  indication  of  how  a  strake  works. 

3.8  Aircraft  model  with  strakes  mounted 
on  its  boattail 


3.8.1  C H  te  r  i^a  f  o r_p  l^ac  i  ng_the_af  ter  body _strakes 

Although  it  will  be  discussed  in  detail  in  the  con¬ 
clusion  section,  integrals  of  DCFKE  for  all  strake  configura¬ 
tions  tested  indicate  that  CFKE  is  lowest  for  strakes  of 

H  =  0.5  and  for  «  =  15°. 

r 

In  order  to  place  strakes  along  the  boattail  of  the 
aircraft,  the  location  of  the  wake  and  its  distance  from  the 
fuselage  had  to  be  found. 

Figure  56a  shows  the  wake  path  along  the  boattail 
area  as  it  moves  downstream.  With  that  information  at  hand, 
a  pair  of  strakes  of  the  same  height/length  ratio  as  the  ones 
tested  were  built  whose  height  when  positioned  on  the  boattail, 
reached  the  afterbody  wake  core  (Fig.  56b). 

The  strakes  were  placed  at  an  a  =  15  with  respect  to 
the  wake  core  and  supported  to  the  nearest  pressure  tap  hole. 
This  allowed  us  to  check  the  pressure  measured  by  the  strake 
and  compare  that  to  the  pressures  obtained  in  the  delta  wing- 
strake  combination  set  up.  Besides,  the  static  pressures  along 
the  boattail  including  strakes  were  measured  as  shown  by 
figure  57. 


Although  the  best  overall  configuration  (lowest  CFKE, 
highest  lift)  would  be  H  =  0,5,  <*r  =  0  strake,  =  15°  were 
used  in  order  to  validate  the  CFKE  criteria  on  which  our 
investigation  was  based. 

The  location  along  the  fuselage  where  strakes  were 
placed  was  arbitrarily  chosen  along  the  region  where  the  pri¬ 
mary  separation  line  was  found  from  the  limiting  streamlines 
vi sual i zation . 

Figures  58a, b  show  how  strakes  were  placed  along  the 
boattail  and  their  angle  with  respect  to  the  afterbody  lines. 

3 . 9  Limiting  streamlines  visualization 

Limiting  streamlines  flow  visualization  of  the  aircraft 
with  strakes  is  shown  in  figure  59.  Figure  59a  shows  how  the 
strake  delays  the  primary  separation  line;  figure  59b  shows 
that  the  downwash  angle  along  the  boattail  is  reduced,  stream¬ 
lining  the  flow,  probably  due  to  the  pressure  field  induced  by 
the  strake.  Figure  59c  compared  to  limiting  streamlines  without 
strakes  shows  that  the  flow  is  attached  to  the  boattail  as  it 
reaches  the  end  of  the  afterbody  area,  eliminating  the  pressure 
loss  and  wake  induced  by  the  afterbody  alone. 

3.9.1  Pressure  measurements  in  the  aircraft  strake 

The  pressure  measured  by  the  strakes  placed  on  the 
aircraft  (Fig.  54c)  agree  with  the  results  obtained  from  the 
pressure  measurements  study;  that  is,  the  strake  feels  a  nega 
tive  lift  on  its  surface.  Its  value  agree  closely  with  the 
ones  measured  on  the  f 1 at-pl ate-del ta-wi ng-strake  set  up. 
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ANALYSIS  OF  RESULTS 


4 . 1  On  strakes 

Besides  the  observations  noted  earlier  on  the  flow 
fields  some  analysis  is  considered  necessary  by  looking  at 
figure  60.  It  shows  the  values  of  the  cross  flow  kinetic  energy 
obtained  for  different  strake  ratios  H  and  different  angles  of 
attack.  For  all  ratios  of  H,  the  CFKE  is  lowest  at  positive 
angle  of  attack. 


The  values  of  CFKE  obtained  with  a  thick  strake  of 
H  =  0.5  are  close  to  those  obtained  with  a  thin  strake.  Also 
the  resultant  flow  field  induced  by  a  thick  strake  is  similar 
in  trends,  and  close  in  magnitude  of  velocity  vectors  to  the 
one  obtained  by  a  thin  strake.  This  allows  us  to  conclude 
that  the  pressure  distribution  on  strakes  and  their  behaviour 
along  a  vortex  is  similarly  independent  of  its  thick  ones. 

4.2  On  the  pressure  distribution  alonq  the  boattail 


The  pressure  distribution  along  the  aircraft  boattail 
with  and  without  strakes  is  shown  in  figure  61.  It  shows  that 
by  properly  placing  strakes  on  it,  the  pressure  on  the  afterbody 
increases.  In  other  words  an  extra  "push"  or  "thrust"  is 
created  due  to  the  effect  induced  by  the  strakes  on  the  surface. 
The  extra  axial  force  obtained  when  strakes  are  placed  on  the 
surface  was  found  by  multiplying  the  wedge  pressure  of  each 
tap  by  the  projected  boattail  area  of  each  hole,  which  was 
chosen  by  drawing  centerlines  along  mid  distance  points  between 
taps.  It  was  found  to  be  23%  extra.  This  force  is  subtracted 
from  the  total  drag  of  the  aircraft.  An  estimated  reduction  in 
drag  for  a  typical  cargo  aircraft  would  be  of  approximately 
0.004  Cx  under  the  conditions  described  in  this  section. 
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Measurement  errors  as  a  result  of  instrumentation 
inaccuracies  as  well  as  test  conditions  fluctuations  were 
deduced  from  calibration  curves,  and  by  monitoring  deviations 
that  occurred  on  certain  test  parameters. 

5 . 1  Accuracy  range  of  equipment  used 

All  pressure  measurements  were  calibrated  against  a 
Betz  manometer  whose  accuracy  is  within  ±  0.1  mm  H2Q.  The 
deviation  of  the  pressure  transducer  used  for  the  probe  and 
surface  pressures  is  of  0.5%. 

Variations  in  the  tunnel < vel oci ty  were  kept  to 
±  0.2  mm  H20  with  set  test  conditions. 

Pressure  variations  in  the  five  hole  probe  tubes  when 
being  calibrated  were  of  about  1.0  mm  H20. 

Angle  settings  in  the  calibration  duct  could  be  easily 
set  within  0.2°. 

5 . 2  Error  margins  in  test  data 

From  the  five  hole  probe  calibration  data  we  can 
deduce  its  roll  and  pitch  angle  deviation.  The  accuracy  of  the 
roll  angle  is  within  2.5°  at  low  values  of  e  and  of  0.89°  at 
larger  settings.  The  probe  misalignment  angle  is  as  low  as  0.4® 
and  a  maximum  of  2.5°. 

From  flow  quality  measurements,  the  flow  angularity 
of  the  wind  tunnel  was  found  to  be  of  about  1°. 


Static  pressure  coefficients  measured  on  the  aircraft 

boattail  and  on  the  strakes  surfaces  varied  within  ±  0.05  of  C  . 

P 


6.  CONCLUSIONS 


The  main  conclusions  which  can  be  drawn  from  this 
study  are  : 

Aircraft  which  have  a  highly  upswept  afterbody  induce 
a  wake-like  flow  which  leaves  the  boattail  as  it  moves  upward. 

Strakes,  when  placed  along  the  path  of  a  vortex, 
streamline  the  flow,  thus  its  cross  flow  kinetic  energy  is 
reduced . 

When  strakes  are  placed  at  positive  incidence,  the 
resultant  flow  field  downstream  of  it  includes  two  vortices 
of  opposite  direction.  By  that,  the  cross  flow  velocity 
components  integral  is  reduced. 

Strakes  provide  lift  when  placed  at  angles  of  attack 
of  zero  or  lower  than  zero  with  respect  to  the  vortex  path. 
Therefore,  a  thrust  is  created. 

When  strakes  are  placed  at  angles  of  attack  equal  or 
greater  than  zero,  a  reduction  in  the  cross  flow  kinetic 
energy  occurs.  Therefore,  the  vortex  drag  is  aiso  reduced. 

Strakes  placed  properly  on  an  airplane  boattail 
reduce  its  overall  drag  by  increasing  the  static  pressure  in 
the  afterbody  surface. 
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TABLE  II 


T1  -  Flow  field  velocity  vectors  tabulation  for  figure  32a 
U  =  2.5) 


T2  -  Flow  field  velocity  vectors  tabulation  for  figure  27a 
(F.S.  6) 


Note  in  all  cases,  figures  are  shown  looking  upstream. 
The  y-axis  increases  to  the  left,  and  z-axis  increases 

upwards . 


The  data  columns  are  tabulated  in  y,z,v,w  order 


O.  281  246E-38  0  0000O0E+O0 

0.0000 0  0  E  +  0  0  O.OOOOOOE+O 0 
0 . 5  0  0  0  0  0  E  +  0 1  O.OOOOOOE+OO 
0  .  1  o  0  0  0  0  E  +  0  2  O.O0OOOOE+OO 
0  .  1  5  0  O  *)  0  E  ♦  0  2  0 . 0  0  0  0  0  0  E  +  0  0 

O  2 0 OOOOE + 02  0 . OOOOO0E+O0 

0 . 250000E+02  0 . 0O0OO0E+O0 

0 . 3 O 0  0 0  0  E  +  0  2  O . 000000E+00 
0 . 3  5  0  0  0  0  E  +  0  2  O . OOOOOOE+OO 
O  .  4  0  0  o  (>  0  E  +  0  2  0  .O0OO00E  +  00 

0 . 45  0  v  0  O  E  +  0  2  0  0  0  0 o  o  O  E  +  0  0 

0 . 5  O  0  0  O  0 £  +  0  2  0.000 O 0 o  E  +  0  0 

0. 5500 0 OE+02  O.000 O00 £+00 

0 . 6O0O0OE+O2  0 . 000OO0E+O0 
0 . 6D00O0E+O2  0.5O0OO0 £+01 

0 . 55  00OOE  +  02  0 . 5OO000E+  0 1 

0 . 5 0  0 0 0 0 £  +  0 2  O.SOOOOOE+OI 
0. 4500 0O £  +  02  0 . 5OO000E+  0  1 

0. 4O0O00E+02  O.SOOOOOE+Oi 
0. 35O000E  +  O2  0 . 5O0O00E+0  1 
0 . 30  0  00  OE  +  02  O . 50OOOOE+  0  1 
0 . 250000E+02  O.SOOOOOE+OI 
0. 2  O 0  0  0  0  E  +  0  2  O . 50 OOOOE+0  1 
0.1500O0E  +  02  <>  .  50O000E  +  O  1 

0 . 1 0  O00O E  +  02  O.SOOOOOE+OI 
O. 5  0  O  O  0  O  E  +  O 1  0 . 50OOOOE+0  1 

0.00O0O0E+00  O . 5O0O00E+0 1 
0 . 00  0  0 O 0  E  +  0 O  O .  1 O0O00E+  02 
O.SOOOOOE+OI  0 . 1 0O000E+O2 
0 . 1 0 O >5 0 O E  +  O 2  0 . 10OOO OE+02 

0 . I5O000E+02  0 . 100000E+O2 

0 . 2  0  0  0  0  0  E  +  0  2  0.100000E+02 

0 . 2500O OE+02  O . 1 00000E+02 
0 . 3  0  O  0  0  0  E  +  0  2  O .  1O0OO 0E+O2 
O . 35000 OE+02  O.100O0OE+02 

0.400O0OE+02  O  10OO0OE+02 
0.45O00OE+02  0 . 1OOOOOE+02 

0.5000O0E+O2  O. 1000O0E+02 
0.5 5 O0O OE+02  0 . 10OOO0E+O2 

0 . 60000OE+02  O . 10OO00E+O2 
0 . 6OOO00E+O2  O  150000E+O2 
O.  55O000E  +  02  0 . 150O0 C'E  +  02 

O. 5OOOO0E+O2  0 . 15O000E+02 
0. 45  0O0OE  +  O2  0 . 1 50000E+  02 
0 . 4  0  O  0  O  0  £  +  0  2  0  .  1  5  O  O  O  O  E  +  0  2 

O . 3  5  O o  o o  £  +  O  2  0.15  O 0 O O  E  +  O 2 

0.3O0000E  +  02  0 . 1  50  000 E  +  0  2 

O.  2 5 O 0 O 0 E  +  0 2  0  .  1  5 O O 0 0 E  +  O 2 

0 . 2  0  0  0  0  0  E  +  O  2  O  .  1  5  0  O  0  0  E  +  0  2 
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0.1I303SE+01  O.638079E+00 

0  .  S88667E +00  -0.148U4E  +  00 
0  .93272  1  6  +00  O . 1 v  2  1 0  2  E  +  0 1 
O  1741996+01  0 . 435435E+0O 

0 . 153439E+01  0. 1 15228E+0O 

0  .  S&77  5  9E  *•  >0  -■>  5423756-01 
0  386.51  56+00  -o  3  1  4  45  1  6  +  01 
0.2342  1  76  +0  1  -  0 . 5 07 Sfc 3 E  +  0 1 
0 . 37  35  3  96  +  0  1  -0. 450  163E  +  01 
Q . 4923476+0 1  - >  1230676*01 
O .  3  6  9  6  5  7  E  +  0 1  0 .211 34  SE  +  01 

O . 21 71 95E+0 1  0.264413E+OI 

0. 1953136+01  0.323937E+01 

0 . 140075E+0 1  0 . 3 15853E+01 

0.148435E+O1  -O . 539769E+00 
0 . 6  3  7  ?  0  3  E  +  0  0  -O .  1  79  41 3E  +  01 
0 .  1 4  3  2  8  3  E  *  0 1  -0.189941E  +  01 

-  O  4  3  8  7  0  3 1  +  0  v  -  O  2  0  5  2  7  8  E  +  0 1 

O.47&261E+0O  V.816345E+0O 

0 .  1 8  430  3t  +  0 1  0 . 2  38  13  3E  +  01 

0.2415746+01  0 . 1 81 182E+01 

O . 2  fa  8  5  0  5  E  +  0 1  0.1O5296E  +  OO 

O. 26 930 4 E +0  1  -  0 .  1  1 006 4 E  + 0 1 
0  .  2  2  8  3  1 5E  +0 1  -O.l 72  31  4E+01 
0.  13  02426  + 01  -V.237049E  +  01 
0.23O55OE-OI  -0 . 252723E+01 

-0 . 36  1935E+00  -  0 .  1 8652 1 E  +  01 
O  0  o  v  o  0  06  +00  — o  .  1  v  4  3  6  1  E  +  0  1 
0  . 3 9  5  9  2  7  6+00  -  0  6  7  0  0  3  fa  E  +  0  0 
0  4  5 fa 5  4  y 6  +<:•  v  -  v  .  6  392706+00 
0 . 4  0  fa  2  5  5  E  +  0  0  0. 580571E-01 

O’  .  fa 5  *  7  3  76  -fa  1  -  0  .  7  5 5  fa  7  2  E  +  0  0 
0  1  5  7  2  4  s  t  *  v  0  -  O’  4  0  8  fa  5  3  E  +  0  0 

-  fa  1  fa  5  s  r  5  c  *  0  -  1  2  7  2  7  7  6  +  0  1 

0  .'252126*  0  >;•  -  0  .  1  3  2  2  8  fa  E  +  v  1 

o  v  ?  0  >.  J  t  + '  ■  o  -O.l  33  56  1  £+01 
Oil 64S26 +0  i  -0 . 95284  1  E  +  00 
0.1590396+01  -O . I 26954E+0Q 
0  .  12  8  8  2  16  +0  1  -  0  4  4  4  2 1  fa  E -  0 1 
v  .  fa  fa  ?  o  ?£  +  0  v  v  .  2  4  7  0 o  3  £  +  0  0 
0  .  fa.' 9  *  fa  v c  - 2  0  2  C’ 7  fa 4  3  t  +  0‘  0 

-  0  .•■  fa  7  7  9  1  c  +  0  0  -  v  3  35  0  5  7  £  +  0  0 

- 0  .  b fa  i  2  2  fa E  +  0  0  -  0  .  5  74  16  ? E  +  00 

0.1816346+01  -0. 493147E+00 


Cassette  digital  data  readout  and  implementation  of 
pressure  transducer  coefficients 

Program  for  operating  data  blocks 

Program  used  to  interchange  data  points 

Program  for  altering  data 

Calculation  of  velocity  vectors,  preparation  of  stream¬ 
lines,  data  matrix,  and  creation  of  pressure  isolines 
block 

Calculation  of  velocity  vectors  and  creation  of  cross 
flow  kinetic  energy  isolines  block 

PET  Commodore  program  used  for  calibration  of  five 
hole  probe 

PET  Commodore  program  used  for  data  acquisition  into  d' 
digital  cassette  (6  and  H  require  machine  language 
program  to  relate  scanivalve  A  at  set  intervals  of 


READ  file  created  by  cox  and  rewrite  it  after  applying  coefficient 


NAME 

AUTHOR 

DATE 


•IE  US 
OTTE  A 
10 -NAY- 8 3 


SUBROUTINES  CALLED  :  PARANI,  GETCDX 


PROGRAM  J  EU S 

DIMENSION  COEFC2  >,  FC<  5  > 
LOGICAL* 1  YES 

NCAN  =  5 

CALL  PAP  AMI  C COEF, 2,1) 

CALL  GETCDX  C NCAN, COEF, FC . FAC T 


SUBROUTINE  GETCDX  C NPTS, COEF  .  PC  .FACT  '• 

DIMENSION  I  VAL<  50  )  .  I B UF  l  <  2 5b  >  ,  I  8U  F  2C  25  6  > ,  C 0  E  r  >  2  >  ,  F 
LOGICAL* 1  F ILE 1C  15  >.F  JLE2C  15  > , YES 
I  V 1  =  1 
I  V2  =  I 

F  AC  T  =  1  .  !  J  EUSETTE 

TYPE  1 

FORMATS'  RAH  DATA  FILENAME  :  '*> 

CALL  GETSTR  (5. FILE!, 14  > 

CALL  ASSIGN  C  1  .FILE1  > 

DEFINE  FILE  1  CO, 256, U,  I  V  1  ) 

TYPE  500 

FORM AT  <  '  NEW  FILENAME  :  '$) 

CALL  GETSTR  C3.FILE2.14) 

CALL  ASSIGN  C2.FILE2) 

DEFINE  FILE  2  < 0 , 256 , U , I V2 ) 


I0UT=O 
TYPE  1111 

FORMATC'  OUTPUT  VERIFICATION  ?  £Y/NJ  :  '*> 

I FC YES <  )  )  I OU T  =  1 

READ  (.  1 '  IV1  >  IBUF1 
WRITE  <  2 '  I V  2 )  IBUF1 

NCAN=18UF 1C  2  > 

I FC NPTS  . GT . 0 )  GO  TO  5 
GO  TO  6 

IF C NCAN. HE. NPTS)  NCAN*NPTS 
NREC  =  IBUFU  3  > 

I FC NPTS . GT .  IBUF1 (2  )  )  NR E C  =  NR EC/NC AN* I B UF 1 < 5  > 

IREC-0 
J  1  =  1 
J2* NCAN 

READ  < 1 ' IV1 >  IBUF1 
I  *0 

DO  25  J-J  I# J2 

I-I ♦  I  - 

EQ.O)  GO  TO  4000  f  JEUS 


35 


COH  T  INUE 
I  RE  C  =  1  RE  C  +  1 
I  GO  =  1 


IFC  IREC.NE.NREO  GO  TO  1000 
IFC IREC . E Q . NREC  >  GO  TO  40 
J 1  =  J2  +  1 

I  FC  </  1  .  EG  .  25?  >  WRITE  C  2  ‘  I  V  2  > 
IFC  J 1 . EG .257)  GO  TO  10 
J2=J 1+NCAN-  1 
IF(J2.LE.256>  GO  TO  15 
1=0 

DO  35  0=J1,256 
1=1  +  1 

I  VALC  I  >= I BUF 1 C J  ) 

CONTINUE 
I  GO  =  2 

GO  TO  1000 

WRITE  C  2 '  I V  2  >  IBUF2 

J  1  =  1 

J2=J2-256 
GO  TO  II 


I  BUF 


K  =  0 
J  E  =  J  2 

IFC  IGO  EGS  .2  >  JE  =  25b 
DO  3000  JJ  =  .J1,JE 
1C  =  K  +  1 

FCC  K  >  =  COEFC 1 > 

DO  2000  11=2.2 

FCC  K  >  =  FCCK.  )+COEF<  I  I  >*FLOATC I VALC K  >  >**<  I  I  -  1  ? 

CONTINUE 

IBUF2C  J  J  )=FC<  K  >/F  A  C  T 
CONTINUE 

IFCIOUT-EO.I)  TYPE  99  9  9  ,  .J  1  ,  J  2  ,  C  I V  ft  L  C  N  >  ,  N=  1 , 5  > .  C  I  B  UF  2(  N  > ,  N»  J  1 ,  J  E  ) 
F0RHA7<2I*.  *  ',516,'  •  '  ,  5  I  6  > 

GO  TO  C60 ,  70  > ,  IGO 

DO  5  000  JJ.J=J.256 
IBUF2C  JJ  J  )  =  0 
CONT  INUE 

WRITE  <  2 '  I V  2  >  IBUF2 
CALL  CLOSE  < 1 > 

CALL  CLOSE  <2  ) 

RETURN 


PROGRAM  APPEND  1 

LOGICAL  *1  NAME  1?  15  >,NAHE2?  15  > • YES 
DIMENSION  IBUFC256). IVALC5> 


TYPE  10 

FORMAT? '  INPUT  FILENAME  t  ' $ > 

CALL  GETSTRC5,NAME1,14> 

DEFINE  FILE  1  <0,256,0..  I V 1  > 

CALL  ASSIGN  C  1  ,NAHED 
I¥l  =  t 

TYPE  15 

FORMAT? '  OUTPUT  FILENAME  ?  '*> 

CALL  GET  STR?  5  .•  NAME  2 ,14) 

CALL  ASSIGN  ?2,NANE2> 

DEFINE  FILE  2  <0,5, U, IV2  > 

I  V2=  1 

CALL  GETVALC5,  IVAL,  ISTAT  >HAHE  > I  V  1  > 

TYPE  1234  -  IVAL 
F0RHAT<5I 10  > 

W  R I  T  E  ?  2 '  I  V 2  > I  V  AL 

IF?  IVAL?  1  >.E(1 ,0>  GO  TO  30 

GO  TO  20 

CALL  CLOSE  Cl) 

TYPE  50 

FORMAT?'  DATA  FILENAME  TO  APPEND  [Y.'Nl 
IF?  NOT  YES?  >  )  GO  TO  55 
TYPE  60 

FORMAT?'  INPUT  FILENAME  :  '$) 

CALL  GETSTR  ?5,NANEI, 14  > 

CALL  ASSIGN  (1, NAMED 
DEFINE  FILE  1  < 0 , 2 5 6 , U , I  V 1  ) 

I  V  1  =  l 

GO  TO  20 

CONTINUE 

DO  56  1=1,5 

IVAL? I  >  =  0 

WRITE?  2'  IV2  >  IVAL 

CALL  CLOSE  <2> 


•  LECTURE  TRAIN  DE  CAHAUX 


TYPE  40 

FORMAT?'  DO  YOU  WANT  TO  CONTINUE  ?  CY/N]  !  '$> 

IF?  YES?  )  )  GO  TO  5 

END 

SUBROUTINE  GETVAL  C NPTS , I V AL , K , F I LE ,  I V 1  ) 
DIMENSION  IVAL? 1  ),  IBUF?  256  > 

DATA  I DEB/O/ 

LOGICAL*!  F  ILE? 1 5  ) 


K*0 

IF? IDEB . NE  0  >  GO  TO  30 
I  V 1  -  1 

INF  I LE  =  0 
DO  2  I«1 ,  15 

IF? FILE? I > . NE . '  ' )  GO  TO  3 

CONTINUE 

GO  TO  4 

INFILE*! 


»  STATUS  INITIALIZATION 


!  FILE  NOT  DEFINED  BEFORE  CAL 


<  NO  BLANC  STRING,  FILE  EXIST 


I  FILE  DEFINED  BEFOR 


mil*!' 


TYPE  1 

FORM ATC / '  RAW  DATA  FILENAME  :  '*> 

CALL  GETS  TR  C  5. F  ILE, 14  > 

CALL  ASSIGN  C  1  .FILE  > 

DEFINE  FILE  1  < 0 . 2 5  * , U ,  I  V 1  > 

READ  C  1  '  I  VI  >  I  8 li F 
NCAN=I BUF  C  2  ) 

I  Ft;  NPTS  .  GT  .  <•  >  GO  TO  5 
GO  TO  '6 

IFC NCAN . NE. NPTS)  N  C  AN  =  N  P  TS 
NREC=IBUFC3  ) 

IFCNPTS.GT.1BUFC2))  NREC =NRE C /NC A N* I 5U F C 5  ) 
I  R  E  C  =  0 
J  1  =  1 
J  2  =  N  C  A  N 

READ  C  1 '  IV1  )  I8UF 
I  =  0 

DO  25  J  =  J  1  /  J  2 
1=1+1 

I VAL< I  ) *  I B 0 F C J  > 

I FC  I VALC  I  )  . ES  . 0  >  GO  TO  40 
CONT I  HUE 
I  R  E  C  =  I  R  E  C  + 1 
I DE  B  = 1 

IFC  IREC.  NE.NREO  RETURN 
GO  TO  40 
J 1  =  J  2+  1 

I FCJI. EC. 25 7)  GO  TO  10 
J2=J 1+NCAN-l 
IFCJ2  LE.256)  GO  TO  15 
1=0 

DO  35  J s  J  I  / 2 5 6 
1=1  +  1 

I VALC  I  >=IBUFC  J  > 

CONTINUE 

READ  C  1 '  I  VI  >  I  BUF 
J  1=  1 

•J2=. 12-25  6 
GO  TO  20 
I  D  £  B  =  0 
CONTINUE 

IFC  INFILE  EQ. 0  >  GALL  CLOSE  Cl) 

K  =  -  1 
R  E  T  U  R  N 


PROGRAM  UPDATE  '  w  *  j -  38  - 

Did  EM  SIGN  IVrtL<5  >,  J  VAL<  3  >' 

LOGICAL*!  NAME1C 15  >.NA«E2C 15  >,YES 
TYPE  20 

F  0  R  H  A  T  <  ■’  NAME  OF  FILE  TO  BE  UPDATED  s  '  $  ) 

CALL  GETS  TP  (5,NAME1<  14  > 

CALL  ASSIGN  <  1  ,  NAME  1  > 

DEFINE  FILE  1  <0,5,0,  I V 1  ) 

TYPE  40 

FORMATS  1  RECORD  NUMBER  TO  BE  REPLACED  s  '$> 
ACCEPT  5  0,  I R 1 
FORMATS  1 5  > 

IF<  IRl  EG  .0  >  GO  TO  999 
TYPE  70 

FORM  AT  < '  NAME  OF  FILE  TO  PICK  UP  RECORD  :  '$) 
CALL  GETSTR  <5,NAHE2,14> 

CALL  ASSIGN  <2,NAHE2> 

DEFINE  FILE  2  <0 , 5 , U,  IV2  > 

TYPE  80 

FORMATS  RECORD  NUMBER  TO  PICK  UP  :  '$> 

ACCEPT  50.IR2 
I F<  IR2 . EG  0  >  GO  TO  883 
REAO(l'IRl)  IVAL 
TYPE  1000, IVAL 
FQRHAT(5I  10  > 

READ<2'IR2)  JVAL 
TYPE  1000, JVAL 
TYPE  2000 

F  0  R  M  A  T  ■'  •'  RECORD  OK?  CY/N3  :  '$> 

I F<  . NOT. YES<  )  >  GO  TO  75 
WRITE*.  1  •  IRl  )  JVAL 
CALL  CLOSE  <2) 

GO  TO  30 
CONTINUE 
CALL  CLOSE  <  2  > 

CONTINUE 
CALL  CLOSE  <  1  > 

TYPE  4000 

FORMATS  DO  YOU  WANT  TO  DELETE  RECORDS  ?  CY/H1 
IF(  .NOT. YES <  )>  GO  TO  5000 

TYPE  4010 

FORMAT< '  NUMBER  OF  RECORDS  TO  BE  DELETED  :  '$) 

ACCEPT  5 0 , N R D E L 
TYPE  4020 

FORMATS  START  DELETING  FROM  RECORD  NUMBER  5  ' 
ACCEPT  50, ISTART 

Jfi=I$TART+NRDEL 
J W -ISTART 

R E A C> <  1  '  J  R  ,  £  H  D  =  5  0  0  0  >  E R R  =  5 0 0 0  >  IVAL 
TYPE  GOOD , J  W , J  R , IVAL 
F  0  R  M  A  T  <  7  I  8  ) 

WRITECl'JW)  IVAL 
J  R.  =  J  R  +  1 
JW*JW+ I 

I F<  I VAL< 1  > . NE . 0  >  GO  TO  4030 
CALL  CLOSE  <  l  > 

TYPE  3000 

F  0  P.  M  A  T  <  '  DO  YOU  WAHT  TO  CONTINUE  ?  IY/H3  s  •$) 
I  F(  YES''  )  >  GO  TO  10 


PROGRAM  BIN  \ 

AUTHOR :  SUS  ORDONEZ  AR '82-83 
PGM  REPLACES  VALUE  FOUND  IN  PAV  DATA  FRO 
TO  VALUE  ASSIGNED  BY  THE  USER 
DIMENSION  J  BUF?  256  > 

LOG  1  CAL  NAMEK  15  ) 

TYPE  20 

FORMAT?'  ItlFUT  FILENAME!'*) 

CALL  GETSTRC5,  NAME  1  ,14) 

CALL  ASSIGN? 1 , NAMEt  > 

DEFINE  FILE  1? 0 , 25 6 , U , I V  1  ) 

I  VI  =  1 
TYPE  3  0 

FORMAT?'  BLOCK  NUMBER  TO  BE  MODIFIED:'*) 
ACCEPT  50, N 
FORMAT?  1  5 ) 

N  *  N  +  1 
NR  *  N 

READCl'N)  IBUF 
TYPE  4  0 

FORMATf'  NUMBER  OF  VALUE  IN  BLOCK:'*) 
ACCENT  5  O , K 
TYPE  200 

FORMAT  C '  ENTER  NEW  VALUE1'*) 

ACCEPT  5  0,  I  V A L 
I  B  0  F  <  K  >  =  I  V  A  L 
WRITE  ?  1  ’ NR  )  I BUF 
CALL  CLOSE? I ) 


FROM  DIGITAL  CASSETE 


sir 


PREPARE  A  DATA  FILE  FOR  S* RE  AML  INES  PLOTTING  PROGRAM  AMD  IS 
N  A II  E  :  NFVSUS 

GUS  ORDONEZ  A  R  '02-33 
DATE  :  2  £  - M  A  Y - ?  ? 

P  G II  CALCULATES  CROSS  F  L  0  M  V  E  L  0  C  I  T  Y  l-  0  '1  P  0  MEAT  S  A  N  D  TOTAL  PR  E 
AND  PREPARES  DATA  FOR  STREmML  I*ES  PLOTTING  AND  iSOLINES  OF 
TOTAL  PRESSURE 
PROGRAM  Nr V GUS 

CO  II II  OH  X  F  \  5  0  )  ,  Vi.  5  v  >  .  tf<  5  0  >  .  I  V  A  L  ■  5  >  ,  PC  5  > ..  XF  L  T  S  0-'.  5  0  > 

COMMON  X  P  R !  5 «  >  ,  V R  <  5  ■>  >  .  w  R  <  5  O  ;•  .  T  f  T  L  E  •'  2  '>  '  HR  -NY 
LOGICAL’''!  N  A  H  E  I  <  1  5  >  <  ii  A  M  E  i  f  i  5  >  NAME  ?•'  t  5  >  ,  V  £  S 


OL  IMES 


S  S  0  R  E 


G  =  9  .  8  0  to  3  5  & 

TYPE  100 

FORMATC  •’  NUMBER  OF  POINTS  OH  X-AXIS 
A CC  E P T  J10, NX 
FORMAT  Y  I  5  > 

TYPE  120 

F  0  R  N  A  T  <  '  NUMBER  OF  POINTS  ON  Y-AXIS 
ACCEPT  110, HY 
TYPE  130 

FORMATY'  INITIAL  VALUE  FOR  X  :  '$> 

ACCEPT  140, XI 
FORMAT.-:  F  t  0  v> 

TYPE  j  5 <> 

FORMATS.  '  INITIAL  VALUE  FOR  Y  :  '*> 

ACCEFT  140, YI 
TYPE  155 

FORMAT  (.  ‘  STEF  FOR  X/Y  :  '*) 

ACCEPT  140, STEP 

TYPE  ISO 

FOR  MATY '  DATA  INPUT  FILENAME  :  '$) 

CALL  GETSTR  Y 5 , NAME  1  >  14  ) 

CALL  ASSIGN  <l,NAMEl> 

DEFINE  FILE  1  <0,5, U ,  I V  1  > 

I  V  1  =  1 


N*MX*NY 
NNM  =  N 
N  I  *  H  +  1 

CALL  ASSIGN  <  2,  'STREAM  OUT 
DEF  INE  FILE  2  Y  0 , 8  ,  IJ  ,  i  V  2  ) 


!  TOTAL  NUMBER  OF  POINTS 


TYPE  500 

FOR  IT  ATi"  '  OUTPUT  FILENAME  FOP  ISOLINES  :  '$> 

CALL  GETSTR  Y  5  ,  NAMES •  14  ) 

CALL  ASSIGN  <  3  •  N  A  M  E  5  > 

TYPE  510 

FORMATY'  TITLE  FOR  ISOLINES  FLOT  :  '$> 

ACCEPT  520, TITLE 
FORMAT ( 20A4 > 

WRITE Y 3, 525)  TITLE  !  RECORD  1  FOR  PLIS 

FOR  MATY  IX , 2  0A4  > 

WRITEY 3,*)  NX , HY  I  RECORD  2 

XPC  1  >*0  . 

DO  530  1=2, NX 
XPY I  >»XP<  I- 1 >  +  ST  EP 

WPITEY3,*)  YXP-.K  ),  K  =  1  ,NX  >  *  RECORD  3  • 

XPY 1  >*0 . 

DO  540  0=2, NY 
XPC  J  >  =  XPC  J- 1 >*STEP 


I  RECORD  2 


*  RECORD  3  -  X-CQORD 


rrrr: 


r:  r.  j  ,  r^r 


n  my  .  -1 


1  oo 
l  •■■  1  <> 


1  0  30 


C  . 


C 

2  0  4  0 


C 


1  8  0 
C 

C 


c 

20 
c  . 
200 

2  1  0 

c 


3  0  0’ 


1111 
4  00 

c 

30 


-  41  - 


E2  * 

•:P.!TE--:  2-  n-  H,0,0  H 

r  y  -r  ■-> 

'*  p  =  y  i  -  s  t  g  p 
T  V  F  t  .1  *  *l  i  v 

F0f?r:  h  T  HO  S  FNEF  I  C  PRESSES  h*  '  $  '• 

h  C  C  E  F  '  14  <;• ,  P  m  T 
TYPE  1020 

F  0  R  M  A  T  C  TEMPERATURE  •  j  :■ 

ACCEPT  1 4  0 ,7 o 
TYPE  l  O  3 

F  0  R  M  h  T  •:  0  Y  HAHIC  PRESS  U  Ft  <  "  m  f-:  2  .3 

ACCEPT  240,&phh 
C  T  E  =  D  P  ff  H  *  G 
XPC  l  -STEP 


D  A  :  ft 


00  30  1=1 .  M Y 

.)  L  -■} 

YP= YP+STEF 
S  I  G  =  1 

IFC  M 0 C- <  JL  .  2  )  .  EG  0)  SIG=-2 
X  M I  M  =  X  ! 

X  PI  A  X  =  3  T  E  P  *  r  L !)  A  r  C  M  X  -  1  >♦  X  I 

I  F  <:  H  0 D  0  J  L  /  2  >  .  £  Q  .  0  ’>  XPC  t  >  =  X  M  ft  X  -'•STEF 

I  F<:  MOD  C  J  L  »  2  >  .  ME  .  0  >  X  P  c  1  >  =  y  M TE  ? 

0  0  2  0  I  =  1  ,  M  X 

XPC  I  >  =  XPC  I  >+STEP*SIG 

IF<  I  .GT.  1  )  XPC  I)=XP<I-1  >  +  S  T  £  F  *  S  i  G 

READ  C  1 '  I  VI  )  I  VAL 

DO  ISO  K  = 1  ■ 5 

PC K  >  =  F  L  0  A  7  C I V  A  L  < K  >  ) 

C  OH  T I  HOE 

CALL  CALCUL  C P , V AR G • W A R G . P A T , TO , B T > 

VC  I  >  = V  AR  G 
WC i )=OARG 
XPRC I >=BT/CTE 

CONTINUE 


DO  210  K=I,NX 

WRITE'  2  '  I  V2  •*  V  PC  K  > ,  YP  ,  we  *  >  ,  w  i  k  v 
CONT [HUE 


I  F  C  H  OD  C  J  L  >  2  ) .  E  0 . 0  )  GO  TO  4<:<v 
K  R  =  0 

DO  300  K  =  NX , 1 , -1 
K  R  =  K  P  *  1 

X  P  L  I  S0(  KR  >  =  X  P  R.  *'  K  > 

WRITE  3<  *  )  <  X  P  L  I  S  0  C  K  >  .•  K  =  1  H  X  >  !  RECORD  5 

TYPE  i  t  1  1 ,  <  XPLIS'K  K  >  .  k>  i  ,  NX  ) 

F0R(1AT<5F1  2  3  ) 

GO  TO  30 

NRITET3/*)  CXPRC  K  ),  IC=t  -  NX  >  !  RECORD  5 

TYPE  1111  c XP Ri  K  >  <  K= 1  . NX  > 

CONT  l HUE 
CALL  CLOSE  Cl) 

CALL  CLOSE  X2  ) . . . . _ : ... _ 


C h lL  CLOSE  <3> 

CALL  ASSIGN  (1,NA«E3> 

CALL  ASSIGN  (2,' ISOLINES  00 

(f :  E  H  D  (  !  •  5  2  5  >  TITLE 

R  E  A  0  <  1  ■  *  >  H  X  .•  H  v 

READ ■  1  -  *  i  t.  \< p  if  >' };'  > ..  »  =  i  4  X 

R  t  A D  v  l  •  *  >  C  XP  R i.  K  > ..  K  = !  .■  N  Y  > 

DO  SOSO  J =1 , NY 
R  E  A  D  <  1  •*  >  <:  XP  L  !S0(  K  >.  K=1  .  NX 
W  R  I  T  E (  2  .•  2  2  2  i  )  XXPLI  SO  <  V  >  •  k  - 
FORM  ATT  1 X , 1  OK  12.5 > 

CONTINUE 
CALL  CLOSE  < 1 ) 

CLOSET  UN ! T  =  2<  D ISP= 'PR INT'  > 

CALL  PEVERS 


DATA  REDUCTION  ROUTINE  FOR  UPDATED  DATA. 

CROSS  FLOW  VELOCITY  COMPONENT'S  AND  TOTAL  PRESSURE  ARE  CALCULATED 
SUBROUTINE  C  A  L  CU  L  <  P  /  V  •  w  /  F  A  7  .  T  *>  >  6  ~  > 

DIMENSION  P<  1  > 

AS  I  NT  X >=ATAN( X/SQRTT 1  . -X*X  >  ) 
ftCOSTX  >=ATAN(  S8RT<  1  ,-X*X  )  /  X  > 

PfIOY  =  <  P<  I  >+P<2  >+P<3  >  +  P(  4  >  >%0  25 
:  (  Z  >-P(.  2  >  +  l  8 
P<  1  >-P< 1  >  +  10 
DD  =  PC2>-P<4> 

EE  *  P<  1 >-P<  3  > 

ADD  =  ft  B  S( DD  > 

AEE  =  A  B  S  (  EE  > 

FF  =  ftEE  -  ADD 

IF  (EE .EQ .0 . >  GO  TO  2 

GG  *  ABS<  DD/EE  > 

FIR  =  ft T AN( G G  > 

FI  *  180 . *F IP  '3.  1 41 59 


FIR 
F  I 

GO  TO 
FI  -  90.0 

FIR  =  3  .  14  159/2 .0 

RA  =  DD*DD+EE*EE 

B I  =  SQR"<  R  A  > 

BIR  =  BI/ASS <PC5> -PM 0Y> 

IF  < FF .LE . 0  .  0  05 >  GO  TO  4 
IF  CEE.GE0.01>  GO  TO  5 

T  1  *  13. 72858*61  P.  +  ( -  2.  112228E-1  >*BIR**2-(  9420 486>*BIR**3+ 

<2  933134E-1 >*81R**4 
T 2  *  -< 2. 69 7529E-2>*BIR**5 

TH  *  T 1  *  T2  -  1  .  24  12  t  7 

D  1  *  (  5  870  1  0 9 E  -  1  >+<4.4&6144E-3  >*tHt(-3.990502E-4>*TH**2 

D2  *  < -5. 3(i2912E-6/*TH**3  +  <  5. 53  8 3  i  3£  -  7  >  *  T  H  •*  + 4 
D 3  *  < -7. 405356E-9>*TH**5 


PD  ■  ( DI  +  D2  +  D3  >*1.0 

Cl  ■  C 7 .065830E-4 >*TH+<-3  © J2390E-4 >*TH**2+( 4 . 834588E-5 > 

*  TH  *  *3 

C2  «  (  -1 . 430393E-6>*TH*«4+<  2. 612108E-8>*TH**3-(2. 319786E-1©  ) 
*TH**6 


10 


c 

ccc 


ccc 

cccc 


8  32 


833 

4 


i 

i 


1 

6 

1 


-  43  - 

E4 

t? H T  *  PAT/ <28?  0*70  :• 

av  =  rtBSC < PC  5 > - P n 0 Y 0*2  . /< PD *RHT > > 

mV  =  SQRTC  AV  > 

8T  =  P<  5  >+CP*PHT*AV**2*0  .*5 

IP  <  <P<  1  >  .GE  .  PC  3  >  >  .  AND .  <.  P(  4  •  .  GE  ?(2)  >f 

IP  <  <  P<  1  )  .GE.  P-:  3  >  >  .  AND.  <  P<  2  >  GE  ft  4  >  f  / 

IF  <<P<1  > .LE. P< 3 > > . AND. < P< 2 > . GE . r < 4 >  >> 
IF  <  <  P  <  1  )  .  L  E  .  P  <  3  >  >  .  AN  D  .  (  F <  4  >  G£.P<2>>> 
TH  =  T  H  *  3 . 14159/ ISO . 0 
IF  C  CP  .  LT  .  0  .  >  CP  =  0.0 
I  F<:  TH  .  LT  .  0.  >  GO  TO  332 
SSTH»SINCTH> 

SSF IR=COS<  FIR  > 

SAL  * .SINCTH >*COS< FIR> 

C  TE  =  1 . -SAL*5AL 

CGft  -COSC  TH  > /SORT  ft C  TE  > 


(?U  =  1  . 
U  U  =  2  . 

eu  =  3  . 
e  u  =  4 . 


CGA  =  cost;  TH  )/<  SQRTC  A8S<  l  . -SAL*SAL  )  >  ) 
CAL=SQRT<  1  .  -SAL* SAL  > 

SGft=$QRT t.  1  .  -CGA* Cuft  > 

TYPE  ,TH,FIR,SSTH,SSFIR,3ftL,CGA,CTE 


FORMATS  TH  .  F  !  -  S  TH  /  CF  I  .•  S  AL  •  C G  A C  T  E  '  /  ?F  1  0 . 4  > 
UU=CAL*CGA  +  4  if 
V-C AL*SG A*A y 

y  -  3  a  i.  *  a  v 

I  F  C  C!  U  EO  .  1  .  .OR.OU  E-3  4  >  V*-y 
IF(  QiJ  E  0  .  I  .  .OR  .  Q'j  .  E  5 . 2  .  )  id  =  -  V 

GO  TO  333 
V  -*0  . 


MaO 
0  U  =  h  V 
RETURN 

IF  (OD.LE  0.01  )  GO  TO  8 

T1  *  2  1  .53483*81  R -13  94*s?*e3  Rw*:? 

T  2  -  +  8  .  3  5  3  0  2  3  *  C  I  F .  *  3  -  3  .2  4  6  7  2  3  *■  8  I «  *  *  4 

T  3  =  (  5 . 3  0  3  7  0  4  E  -  1  > »  £  I  ft.  «  «  5  -  4  .  i  4  8  4  7  0  £  -  2  >  *  B  i  R  *  *  t 

TH  =  T  1  +  T  2  +  T3-  .7471844 

D 1  =  C  3 . 9093289E- 1 >♦<  8 . 1 93470E-4  >*TH-<6 . 7776  13E-5  >*TH**2 

02  =  C -3  96  041 5E-6 J*TH**3+< -1  . 53  78 24E-?  >*TH**4 

0  3  35  i  1  i'5385'fE'S  -  ?  3  4  7  <  8  <>  £  -  1  o  >  *  T  H  *  *  8 

PD  =  (  0  1+02  +  03  >*  1  0 

Cl  =  1 '  -  5  .  2  5  i  4  6  8  E  -  3  >  *  T  H  «- 1  i  2  i  3  6  5  7  E  -  3  >  *  f  H  *  *  2  - f.  1  .  4  3  8  8  3  £  -  4  J  * 

T  H*  *  5 

C  2  =  (8  7S8080E-6  >*TH**4r»  -  2 .  1  7 250 9E -? >*7H**3* <  t . 907005E-9) 

*  T  H  *  fc  £ 

CP  «  C 1  *  C  2  *  9  ?7,586E-4 


GO  TO  jo 

T  1  a  i  7  2  £  £  *  4  *  R  T  p  —  7  5'5£  2-i  ■*£:»»■»  2 

T  2  *  5  .  7  4  8  9  *>  9  *  8  i  R  *  *  3  -  «  2  i  4  8  7  4  2  < »  8  t  ft  *  *  4 

T 3  a  C 3 .965276E-1  >*eiR**3-< 2  3 Z i u v 3E - 2 > *E I R* *6 
TH  =  T  I  +T2+ T3+0 . 494 192 1 

01  =  <  5  .  95200E-1  >♦<  7  .  554778E-5  ;*  TH-m;  i  .564457E-4  >*TH**2 

02  a  C -4 . 345395E-3>*TH**  3-m  2  1  £  2  5  4  0  E  -  8  >  *  T  H  *  *  4 

03  *  ( -4 . 38095  IE-8 >*TH **5+<  3  248  74  IE- 10  >* TH**6 

PD  »  ( D  1  +  D2  +  D3  >* 1  0 

Cl  =  (-2. 83S8€8E-3>*TH*<  6  1  9 9 82 0 E  -  4 > * TH *  *  2 -C  3 . 4 55 3 7 8 E- 3 > * TH ** 3 
C  2  =  (  1 .431172E-4 >*TH**4-i  2  517231E-8>*TH**5+l 7.831602E-11) 

*TH**8 

CP  *  C 1 ♦C2* I  .  4 0 2 8  8 5 E -  3 


GO  TO  10 

T  1  •  12  .  12255*BIR*<2. 329581  >*< 81 R**2 >-3 . 229866*6 IR**3+ 

C  1 . 284068  >*BIR**4 


•  t  _  *  KDK1C.I 


1 a  aaa  »  g . 


TH  *  T 1  +  T  2  -  0 

D  1  =  <5  .  83329SE-1  >  +  C  2 . S77764E-3  >*  TH  +  C  -  1 . 635622E-4  >*TH**2 

D 2  *  (-  1  .  4225  16E-5  >*TH**3  +  <  6  ?2Si?93E- 7  >*TH**4 

0  3  *  C-8  .S530S1E-')  >*TH**5 

PD  =  < D 1 +02+D3 >*1  00 

Cl  *  <  -8 . 130340E-4  >* TH -C  2  .  0&68S3E-5  >*TH**2+<2  829430E-5  >*  T  H*  *3 
C 2  =  <-*  .  975929E-7 >*TH**4+f  2  329037E-9 >*TK**5 

CP  -  C1+C2+7 . 536941E-4 
GO  TO  10 
END 


SUBROUTINE  USED  TO  PLOT  CROSS  FLOW  VELOCITIES  ON  PDP  PLOTTER 
SUBROUTINE  REVERS 

C  OH  HON  XPC50),  VC  50  >  ,  WC5  0  >,  IVALC  5  >,  PC  5  >,  XPLISOC  50  > 

C  OH  HON  XPRC  30  ), VRC  30>, WR<50>, TI TLEC20 ) , NX ,NY 

CALL  ASSIGN  <2, 'STREAM  . OUT'  ) 

DEFINE  FILE  2  <0,8, U, IV2 ) 

IV2*t  - . - 

CALL  ASSIGN  ( 3 , ' ST R EA HP . OU T *  ) 

DEFINE  FILE  3  <0,8,U, IV3  > 

I  V3  3  1 

READC  2 ' I V2 )  N , NZ , RNZ 
WRI TE< 3'  I V3  >  N , NZ, RNZ 

DO  40  J  =  1  , N Y 
DO  10  1*1, NX 

READ(2'IV2>  XP< I  ), Y , V< I  ) - W< 3  > 

CONTINUE 

L*0 

DO  20  K=  NX , 1,-1 
L*L+  1 

VR<  L  )  =  -V<  K > 

WRC  L  >  =  WC  K  ) 

CONT  INUE 

DO  30  X* 1 ,NX 

WRI  TEC  3'  I  V3  >  XPC  K  >  >  Y  >  VRCK  >,  WRC  K  > 

CONTINUE 
CONTINUE 
CALL  CLOSEC  2  > 

CALL  CLOSE < 3) 

RETURN 


O'- t-PARS  A  DATA  FILE  FOR  ISOLIMES 

NAME  :  14  E  U  PROSE  AND  VELOCITIES  KINETIC  <  iiFV !'  I  NET 

MODIFIED  FOR  CROSS  FLOW  KINETIC  ENERGY  INTEGRAL  -  C . CRDOHE 
DATE  ‘  2fc-MAY-8? 

POM  PREPARES  DATA  FOR  KINETIC  ENERGY  ISOLINSi  AND 
CALCULATES  CP  OS  FLOW  KINETIC  IN  '"ORAL  NONOI MENS! 0  HALT  ZED  N 
MAINSTREAM  VELOCITY 
PROGRAM  NPVKINET 

DIMENSION  KPf  50> , V<  50  >  ,  IK  5  0  *  ,  I  VAL  ■  5  I-  ,  P  C  5  >  .  SPLISOC  50  CFKEL 
DIMENSION  XPRC  50  >  VRO  50  )  N  R*'  5  o  '•  •  ~  !  ~L  E*  2  •>  >  .  D  C  Ffc  EL  <  50  > 
LOGICAL*!  NAME U  15  )  , NAME2T  15  >  N AN E 3<  I  5 > , Y ES 


RHOAIR®! .225 
R  H  0  W  A  T  *  1  0  0  0 . 

G*9.8<>«336 
TYPE  100 

FORMAT.;--  NUMBER  OF  POINTS  ON  X-AXIS 
ACCEPT  120, NX 
F  0  R  H  A  T  <  15 ) 

TYPE  120 

FORMATS  HUMBER  OF  POINTS  ON  Y-AXIS 
ACCEPT  110, NY 
TYPE  ISO 

FORMAT  v  •'  INITIAL  VALUE  FOR  X 
ACCEPT  140, X I 
FORMATCFIO.O) 

TYPE  150 

FORMAT.''  INITIAL  VALUE  FOR  Y 
ACCEPT  140, YT 
TYPE  155 

FORMAT.:  '  STEF  FOR  X/Y  »  '$> 

ACCEPT  140, STEF 


'  S  > 


*  > 


TYPE  160 

FORMAT.''  DATA  INPUT  FILENAME 

call  gststr  <:  5 ,  name  i  >  14  > 

CALL  ASSIGN  <  I , NAME  l  > 

DEF  I  HE  F  T  L E  1  >'0.5.  U,  I V  1  i> 

I  V  1*1 


'  %  > 


'  $  ) 


'  *  > 


tf  =H  X  »«Y 
H  !•)  N  ~  M 

N  1  *  H  +  l 


0  ■<  L  ri  J  •*  r,  “  R  i)  1 


T  v P  E  ■■ 

~  *• » A  A  '  •  ••  fspT  F'LEH*»“£  FOP  ISOL  IN?  S  *  *  ? 

C  At.  L  G  E  *  S  T P  1  5  •  N  A K E  3 ,  1  4  > 

CAL-  ASSIGN  «  3, NAME  3 > 

CALL  ASSIGNS.  '  K.  I  HER  .  OUT  '  > 

DEF I  HE  FILE  4  f 0 , 8, U,  IV4  > 

I  V4  =  2 
TYPE  5 1 0 

FORMAT''  TITLE  FOP  ISOL INFS  PLOT  ;  ' * > 

ACCEPT  520, TITLE 
F  0  R  H  A  T  <  2  0  A  4  > 

WRI  TEC  3,  525  )  TITLE  1  RECORD  1  FOR  PL  IS 

PORMATC 1 X , 2 0 A  4  ) 

WRITE  <3. *>  NX,  NY  1  RECORD  2 

X  P  •:  1  >  *  0 
DO  530  I *2, NX 

XPC  I  >»XP<  1-1  >  +  STEP  .... 

WRITE (3,*>  <XPCK),KS1,NX)  !  RECORD  3  -  X-COORD. 

_  ■  -  ■  ■- _ _ _ _ _ _ _  -  -  '  -  l  i  I  in  Hi  —  l 


DO  540  <J  =2 .- NY 
XP*  J  >  =  XP*  3-  1  )  +  STEP 

WRITE*.  3  /  *  )  <XP*K>,K  =  1,HY>  !  RECORD  4  -  Y-COORD. 

YF=Y i-STEP 
TYPE  1010 

FORMAT*'  ATMOSPHERIC  pressure  .  t  > 

ACCEPT  140, PAT 
TYPE  1020 

FORMAT*'  TEMPERATURE  :  '$> 

ACCEPT  1 4  0 , T  6 
TYPE  1030 

FORMAT*'  DYNAMIC  PRESSURE  <HH.H20>  i  $) 

ACCEPT  1 40, DPMH 

VELSQ'C  2  .  ♦RHOWhT*DP«H  > PH  Oh  I R* i 0 0 0 .  > 

CTE«DPf1H*G 
XF*  1  )» XI -STEP 

DCF  RET *0  . 


DO  30  J=1,HY 
■J  L  =  I 

YF*YP+STEF 
S  I G  -  1 

I F( MOD* J L , 2 > . EC . 0 »  S I G- -  1 
X  H  T  M  -  X  I 

XMAX«STEP*FLDATrNX-l  )  +  XI  - 

I  F*  MOD*  tl  L  »  2  )  .  EG  .  0  >  XP< 1  >=XP»hX  +  S T  E  P 
IF*  MOD*  4  L  ,  2  ) .  ME  .  0  )  XP*  1  >«XH3  N-STEr 

DO  2  0  1=1  ,  N  X 

XP*  I  )=XP*  I  >+STEP*S I G 

IF*  I  . GT .  1  >  XP*  I  )  =  X  P  *  I  - 1  >tSTEP*S  £G 

READ  *  1' I  VI  >  I  V  A  L 

DO  ISO  K  =  l ,  5 

PCK  >=f LOAT*  IVALC  K  )  ) 

COST  I  X  U  E 

CALL  CALCUL  * P » VARC . WAR G  .  ? AT  .  Ti  > 

V*  I  )  =  V4«  G 
0*  I  >  *  W  A  P.  G 

CPKEL*  I  >  -  <  Y  •.  I  )  *  *  2  +  W  *  1  vo*  2  > 

D C F  K  c !.  *  T  )  =  C  F  <*  E  L *  I  V<  VEi.  5  0  • 

D C F  K  £  i  -  •>  •'  r  t  "  *  !> 0 r  R  EL*  J  > 

TYPE  lilt-  D  0  K  E  ’ 
n  P  R  <  I  )  =  D  G  F  I,  EL*,  i  > 

CONTI  N  U  E 

I  F  *  M  GO  *  4  L  .  2  >  .  E  Q  .  0  >  GO  T  G  -  o .. 

K  ?  •-  o 

D  G  3  v  v  (*  -  H  X  ,1,-1 
r  R  =  K  R  f  1 

;:pl  iso*  kr  >=:<pp.< k.  > 

WRITE*  3,*  5  *XPl  ISn«  R >,K=1 , MX  i  I  RECORD  5 

00  2  5>  I  =  i  • NX 

WRITE  *  4  ••  IV4  )  V*  1  ) ,  y  ;  I  ) ,  DCF  REL*  i  ) ,  DCF  RET 

CGH  f  I;IUE  ..  . 

TYPE  111!-  *  V  *  I  ),(*■:  I  >,DCFKEL*  I  >»VELSC,DCFKET,  I  « 1  «  N  X  > 


.  .1  s  4 


i  7  »• 

.‘.7  7 


4  ) 


.)  *■  t  i  r  » 


i  i? :  w 


0  ! 
f.  ; 


D  j 
F  l> 

C  1 

T  H  *  *  3 
C  2  * 

t  T  H  *  f  <} 
CP  = 
GO  TO 
T  1  = 


,•  s  c  _  i  t  7 1  h  . 

*  —  ?►••••»  i  ;■  f‘  -  ^  .  ».  n  •  • 

(  >  •  .')!■  •  i  f.  -  •  »  r  f-.  ,  r  \ 

•:  t*  i  +>■•£  +  i>  7  »  »•  i 
1  —  j  !j  i  li  -j  •'  •*  J  ,<  »  T  H  r  ’ 


<  8  7  *  *)  •*  4  *:•  F  -  6  :•  *  r  h  f  * .» 
Ci+C  i*1*  77735bE-4 


. :  :• r  h  - :  ;  4 
7  !  7  if  *  *  T  H  -M  1 


T2 
T3 
TH 
:•  i 
E>2 


17  .  2  ‘>fe  9  4*31  R  -3 . 7 3 30 24  -»3!  R-*  -  : 

5  7  4  £?'■'?  »  5  I  *?  +  *?-  f  2  .  ;  4  r  ••  r  •  * «  ,  »  ,  T  4 

<  3  .9  *?  5  2  7  6  E  -  1  '•  ♦  P  f  P  *  *  'j  -  <  .?  i  •  7  E  -  2  » *  p  !  R  »  *  *j 

T  1  ♦  T  2  ♦  T  3  *  0  4  4  4  1^21 

(  3 .  93200E-I  >•*•»'  7  5  ■>  '#  .’  r  a  £  -  '3  ■  T  H  + 1  l  5*443?E-4>*TH*+2 
<-4.34  5S9,5fc-3>*THi‘*s*»2  l-s2  34.>E-6>MH*  +  4 

<  -  4  38  0  95  l  E  -  3  )  *  T  H  ♦  *  C  3 . 2 4 1>  74  1  E  -  l C  >*  T  H*  +  6 


in  ijj 


V.V.  . 
V.V  ,% 


r>y. 

H*  *  •  *  m 


>'S 
%  s  * 

s'  O 
*.  •*.  A 


kfr  ^ 


C'.C'-C* 
r -.*  v*  . 

*W 

>  .*  .• 
•*  V. 


AV>-* 


kv£ 

V*  1 


j* 


cv:- 

.S’  • 


1  OPEN  1,4  \J\ 

2  OPEN  4,4,4 

5  POKE  348,5 

6  POKE  351,1 

7  POKE  353,1 

8  POKE  852,4 

10  PR I NT #4 : REM  ENABLES  ERROR  DIAGNOSTICS 
20  B I M  I  <  3  ? ,  X  <  5  > ,  V  C  5  ) ,  P  (  5  > 

30  DIM  CF<5>,VC<5> 

55  ST  < 1 >  =0 
30  READ  R 
70  DATA  287 

SO  FOR  1=1 T05 • READ  CFC I >  ■ NEXT 

90  DATA  133. 32, 9. 806356, 9. yu6y5E+u4, 8o94. 3, 1 

100  I NPUT " DATE " I  < 3 > ,  1(4 > ,  1(5 > 

110  POKE  852,4 

120  INPUT "NUMBER  OF  TRANSDUCERS  TO  CALIBRATE, 0  TO  END"; NT 
130  IF  NT =0  THEN  1000 
140  POKE  350,1 

180  I NPUT "WHEN  READV  TYPE  1“;KK 
190  IF  KKOl  THEN  130 
200  HP=0 

210  I NPUT "VALUE  OF  POINT  AND  UNIT  INDEX" ; RP, IN 
220  IF  IN=0  THEN  325 
230  NF-NP+1 

240  REM  CLEAR  AND  INITIALIZE  STORE 
250  SVS ( 320 17) 

280  Y<NP)=RP#CF<IX> 

3 1 0  V=PEEK <  380  >  :+:258+PEEK <  88 1  >  -2048 
320  X<NP>=V -PRINT  V  :  GOTO  210 
325  GOSUE  5000 : REM  LINEAR  REGRESSION 
330  PR I NT# 1, "CALIBRATION  OF  CHANNEL" 

340  PR I NT# 1 

350  PR I NT# 1 , " PRESSURE " , " D I G I TS " , " CALCULATED  PRESSURE " 

355  ER=0 

380  FOR  J-1TO  NP 


363 
365 
370 
380 
400 
410 
420 
440 
450 
1 000 
1012 

1013 

1014 

1015 


VC  <  J  >  = A0+ A 1 #X  < J  > 

ER=ER+ ABS  < VC  <  J ) -V  <  J  > >  * 1 00/V 1 

PR  I  NT#  1 ,  V  <  J  > ;  "  PA  " ,  N  <  J  > ,  VC  <  J  > "  PA  "  :  NEXT 

PR I NT# 1 

IF  ERC1.00  THEN  420 

PR I NT# 1, "POSSIBLE  ERROR  IN  CALIBRATION  OF  CHANNEL 
PR I NT# 1 , " REGRESS I ON  COEFF I C I ENTS  " ; A0 , A 1 
PR I NT# 1 , " MEAN  ERROR " ; ER : PR I NT# 1 
GO  TO  110 
POKE  850,0 

I  NPUT "  P-ATM ,  CONV :  1  FOR  HG,2  FOR  H20".:PAT,IC 
PAT=PAT*CF< IC> 

INPUT"P-DVN  IN  MM  H20";P1 
P1=P1*CF<2> 


II 


1017  I NPUT "ANGLE  OF  ATTACK" ;AL 
1020  INPUT  "TRAVERSE  A";  K2> 

1030  POKE  02,123 

1031  POKE  01,123  : REM  INI  USR  ADD 

1032  INPUT  "TEST  NUMBER"; I<1> 

1033  FOR  J=1T05 

1034  A=USR<I(J>> 

1035  NEXT  J 
1040  A=USR  <  - 1 ) 

1090  I NPUT "TOTAL  TEMP  IN  DEG-C";T0 
1092  T0=TO+273. 15 

1095  RHT=PAT/<R*T0> 

1096  U1 =SGR  <  2$P 1 /RHT  > 


'  J"  IT  /  ^  . 


1200 

1202 

1203 

1204 

1205 

1206 

1207 

1208 


-  51  - 

PR  I  NT#  1 , 11 BATE " .:  I  (  3  ) ;  1(4 ) .:  I  (5  > :  PR  I  NT#  1 
PR I NT# 1 , " TESTS  COND I T I OHS " : PR I NT# 1 
F'RINT#1 ,  "ATM  PRESSURE" ;  PfiT;  "PH" 

PR I NT# 1 , " TEMP  TOT ALE " ; T0 ; " K " 

PR I NT# 1 , " TRAVERSE  A " ; I ( 2 ) ; " FS " 

PR I NT# 1 , " TUNNEL  SPEED " ; U 1 ; " M/S " 

PR  I  NT#  1 ,  "  PR .  DVN  " P 1 ;  "  PA  " 

MU=1 . 4962E-06#T0fl . 5/(T0+120> 


G  2 


1 209  RE=RHT*U 1 * . 062/MU 

1210  PR I NT# 1 , " RE  NO . " ; RE : PR I NT# 1 

1211  PR I NT# 1 , " TEST  NUMBER " J I ( 1 ) : PR I NT# 1 : PR I NT# 1  -  PR I NT# 1 : PR I NT# 1 : PR I NT# 1 
1290  POI=0 

1 300  I  NF'UT  "  EQU I L I  BRAT  I  ON  DELAV  " ;  ED 


1310 

POKE  849, ED 

«.v 

1320 

POKE  59459,2 

1330 

POKE  59471,1 

1345 

POI=F’OI  +  l 

V 

V 

1350 

FOR  J=1T05 

I 

1370 

SVSC32017) 

1380 

POKE  59471,4 

1398 

POKE  59471,1 

cT 

1400 

SVS  C 3 1 945 > 

\v 
1 . 

1430 

V=F‘EEK  (  8bU  )  $ 

148S 


1440  P(J)=A0+A1*V 

1442  A=USR(V> 

1443  NEXT 

1 446  PR  I  NT#  1 ,  "  PO I  NT  " ;  PO I ,  "  XS  " ;  XS ,  "  VS  " ;  VS  •'  PR  I  NT#  1 
1 455  PMOV= < P < 1 >  +  P ( 2  > +P ( 3  > +P < 4  > > /4 
1 460  PDF- ( P ( 5 > -F'MOV > /P 1 
1465  DD=P ( 2 ) -F' ( 4 )  HH=P(1>-P(3) 

1 468  ADD=AES  (  DD )  ■'  AEE=ABS*  (  HH )  •'  FF= AEE- ADD 

1469  IF  DD=0  THEN  1475 

1470  GG=AES<HH/DD>  : F I R= ATN ( GO > :  F I  =  1 88$F I R/3 .14159  ■ GO  TO  1480 
1475  FI =90 

1480  RA=DD*DD+HH*HH 
1484  EI=SQR(RA> 

E I  R=E  I  /AES  (  P  ( 5  >  -F'MOV  > 

CP=1-(P(5>/P1 > 

IF  P  C 1 >  >=P  C  3 ) ANDP ( 4 ) >=P < 2 )  THEN  QU=1 
IF  P ( 1 ) >=F‘ ( 3 > ANDP ( 2 > >=P ( 4 ;• THEN  QiJ=2 
IF  PCI) <=P ( 3 ) ANDP < 2 > >=F ( 4 >  THEN  GU=3 
IF  P ( 1  > OP ( 3 ) ANDP ( 4 ) >=F' ( 2 >  THEN  GU=4 
PRINT#  1  .•  "F'l " ;  PC  1  > ;  "PA" ,  "F'2" ;  P(2> ,  "P3" ;  P<3> ;  "PA" =  PR  I  NT#  1 
PR  I  NT#  1 ,  "  P4  " P  <  4  > ;  "  PA  " ,  "  P5  " ;  P  C  5  > ;  "  PR  " "[FI]".:  FI;  "  CP " CP  =  PR  I  NT#  1 
PR  INT#1,  "EIR"  :BIR,  "  PDF' " PDF’ ,  "  F'MOV  " ;  F'MOV  :  PR  I  NT#  1  •  PR  I  NT#  1  :  PR  I  NT#  1 
GO  TO  1850 
GO  TO  1525 

PR  I  NT  "  ACQU I S I T I  ON  SUSPENDED .  FOR  NEXT  PO  I  NT ,  TVF'E  1 " 

INPUT  KK 

IF  KK=1  THEN  1310 
SVSC31862) 

PR  I  NT  "  ACQU  I S I T I  ON  SUSPENDED.  TO  RESUME,  TVF'E  1" 

INPUT  KK 

IF  KK=1  THEN  110 
GO  TO  2200 
END 


1486 

1527 

1528 

1529 

1530 

1541 

1542 

1543 
1550 
1710 
1850 
I860 
1870 
1900 
1910 
1920 
1930 
I960 


5000 

5010 

5020 

5030 

5040 

5050 

506© 

5070 


XI =0 : VI =0 : XV=0 : X2-0 : V2=0 
FORJ=lTONP 

X 1 =X 1 +X ( J ) : V1=V1+V( J) 

XV=X V+X <  J  )  #  V  <  J  > 

X2=X2+X< J> 12 ■ V2=V2+V< J)  T2  =  NEXT 
DN=X2'-X  1  #X  1  /NP 
A 1  =  <  XV-X 1 # V 1 /HP  > /DN 
R0= ( V 1 -fl 1 #X 1 > /HP 

oo-Qi  *ai*/uQ.mn(mvup.^fw-v  1  *7iyuPj 


s’/, 


\  ' .  • 


.  s.  .  . 
**' '-V- 


:  \  s  -m 


•v  •  . 

V  • ' 


;;/n>i 


:#Si 

r  •  » 

■  Y  >y 


‘iii 


.  .% 
if  : 


•  •'*/ 


!•■  Cl 


•  •  •  •  ■ 


.> 


A 


V 
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5100  RETURN 

8000  I NPUT " NUMBER  CHANNEL " ; N 
8020  POKE  40970,0 
8030  POKE  40971,0 

8040  fl=PEEK< 40972 > : IF  NOTH  HMD  128  THEN 
8050  fl=PEEK< 40973 > : IF  NOTH  AND  128  THEN 
8060  A=A-248 ; GO  T 0  8038 
8078  H=fl+8 

8080  A=A#256+PEEK  4 89 7 4  )  -2048 
8090  PR I NT “ ON  CHANNEL  NUMBER " ; 1 ,  A ; " B I G I T 
8108  GO  TO  8080 
READY . 


1  OPEN  1,4  M 

2  OPEN  4,4,4 

5  POKE  843,5 

6  POKE  851,1 

7  POKE  853,1 

8  POKE  852,4 

16  PRINT#4 : REM  ENABLES  ERROR  DIAGNOSTICS 
20  DIM  I < 8 > , X  < 5  > , V < 5  > , P ( 5  > 

30  DIM  CF ( 5 > , VC < 5 ) 

55  STa>=0 
60  READ  R 
70  DATA  287 

80  FOR  I = 1 T05  •'  READ  CF  <  I  >  •  NEXT 

90  DATA  1 33 . 32 ,  ‘3  •  806356 , 3 . 8O665E+04 , 6834 . 8 , 1 

1 00  I  NF'IJT  "  DATE  " ;  I  < 3  > ,  I  <  4  > ,  I  <  5  > 

110  POKE  852,4 

120  INPUT "NUMBER  OF  TRANSDUCERS  TO  CALIBRATE, O  TO  END", NT 
130  IF  NT =0  THEN  1000 
140  POKE  850,1 

180  INPUT "WHEN  READV  TYPE  1";KK 
130  IF  KKOl  THEN  180 
200  NF-0 

210  INPUT" VALUE  OF  POINT  AND  UNIT  INDEX";RP,IX 
220  IF  IX=0  THEN  325 
230  NP=NP+1 

240  REM  CLEAR  AND  INITIALIZE  STORE 
250  SVS<32017> 

280  V<NP)=RP*CF<IX> 

3 1 0  V=PEEK < 860  > +2 56+ PEEK <  86 1 ) -2048 
320  XCNP)=V ; PRINT  V : GOTO  210 
325  GOSUE  5000 -REM  LINEAR  REGRESSION 
330  PR I NT# 1, "CALIBRATION  OF  CHANNEL" 

340  PR I NT# 1 

350  PR I NT# 1 , " PRESSURE " , " D I G I TS " , " CALCULATED  PRESSURE " 

355  ER=0 

360  FOR  J=1T0  HP 
363  VC < J ) = AO+ A 1 #X < J > 

365  ER=ER+ ABS < VC ( J  > -V < J > > # 1 00/ V 1 

370  PR I NT# 1, V ( J > ; "PA" , X< J> , YC< J) ; "PA" : NEXT 

380  PR I NT# 1 

400  IF  ERC1.00  THEN  420 

410  PR I NT# 1, "POSSIBLE  ERROR  IN  CALIBRATION  OF  CHANNEL" 

420  PR I NT# 1 , " REGRESS I ON  COEFF I C I ENTS  " ; AO , A 1 
440  PR I NT# 1 , " MEAN  ERROR " ; ER  =  PR I NT# 1 
450  GO  TO  110 
1000  POKE  850,0 

1012  INPUT" P-ATM, CONV : 1  FOR  HG,2  FOR  H20" ; PAT, IC 

1013  PAT =PAT#CF < I C  > 

1014  INPUT" P-DVN  IN  MM  H20";P1 

1015  P 1 =P 1 #CF  <  2 ) 

1017  INPUT "ANGLE  OF  ATTACK" ,AL 
1020  INPUT  "TRAVERSE  A"Jl<2> 

1030  POKE  02,123 

1031  POKE  01,123  :  REM  INI  IJSR  ADD 

1032  INPUT  "TEST  NUMBER" ;  1(1) 

1033  FOR  J=1T05 

1034  A=USR<I<J>> 

1035  NEXT  J 
1040  A=USR<-1 > 

1090  INPUT "TOTAL  TEMP  IN  DEG-C";T0 
1092  T0=T0+273. 15 

1095  RHT =PAT / ( R#T0 > 

1096  U1=SQR<2*P1/RHT> 


12UU 

1202 

1203 

1204 

1205 

1206 

1207 

1208 

1209 

1210 
1211 
1290 
1300 
1310 
1320 
1330 
1345 
1350 
1370 
1380 
1390 
1400 
1430 
1440 

1442 

1443 
1446 
1455 

1464 

1465 

1468 

1469 

1470 
1475 
1480 

1484 

1485 
1490 
1495 

1497 

1498 
1510 

1515 

1516 

1517 

1518 

1520 

1521 

1522 

1525 

1526 

1527 

1528 

1529 

1530 

1531 

1541 

1542 

1543 

1544 

1545 
1550 
1600 
1601 
1602 


PR  I  NT#  1..  "DATE" ;  I  C3> I  C4> I  (5) :  PR  I  NT#  1  * 

PR IHT#1," TESTS  CONDITIONS" : PRINTttl  ,  _ 

PRINT# 1, "ATM  PRESSURE"; PAT; "PA"  H  / 

PRINTttl, "TEMP  TOT ALE" ; TO; "K"  1  1  L~ 

PR I NT# 1 , " TRAVERSE  A " ; I C 2  > ; " FS " 

PRINT# 1 , " TUNNEL  SPEED " ; U 1 ; " M/S " 

PR I NT# 1 , " PR . DVN " ; PI ; " PA " 

MU=1 . 4962E-06#T0T1 . 5/CT8+128)  « 

RE=RHT*U  1  * .  062/MU  ! 

PR I NT# 1 , " RE  NO . " ; RE : PR I NT# 1 

PRINTttl , "TEST  NUMBER" ; I C 1 > : PRINTttl : PRINT# 1 : PRINTttl : PRINTttl : PRINTttl  - 
POI-0 

I NPUT " EQU I L I BRAT I ON  DEL A V " ; ED 
POKE  849, ED 

POKE  59459, 252  REM  VIA  DDR  PINS  8,1  INPUT  • 

POKE  59471, 12 : REM  PUTS  1  TO  OUTPUT  PINS  248  : 

POI-POI+l 

FOR  ,T=1T05 
gyc  ( j  *i 

POKE" 59471 , 4 : FOR  K=1 T095 ' NEXT 
POKE  59471 , 12 : FOR  K=1 TO580 = NEXT 
SVS< 31945) 

V=PEEK  (  860  )  *256+F'EEK  C  86 1 )  -2848 

PCJ)*A0+fll#V 

A=USRCV> 

NEXT 

PRINTttl , "POINT"; POI , "XS" ; XS, "VS" ; VS : PRINTttl 
PMGV= C P  C 1 > +P  < 2 ) +P C 3 ) +P  C  4 > ) X4 
PC2)=PC2)+18  : PC  1 >=PC 1 >  +  10 

DD=P  <  2  >  -P  (  4  >  ••  HH=P  C 1  >  -P  C  3 > 

HDD- ABS  C  DD  >  ;  AEE=AESCHH)  :FF=AEE-ADD  • 

IF  HH=8  THEN  1475 

GG=ABS C DD/HH >  : FI R= ATN ( GG > :FI=1 80#F I R/3 .14159  ■ GO  TO  1 480 
FI  =90 

RA=DD*DD+HH*HH 
B I =SQR C RA  > 

B I R=B I X ABS C P  C 5  > -PMOV > 

IF  FF<=0. 805  THEN  1650 
IF  HH>=0. 01  THEN  1600 

T 1 = 1 3 . 72858*3 I R- < 2.11 223E- 1 >  *B I R  T2- C 8 . 942849 >  *B I R  13+  < 2 . 933 1 34E- 1 > *B I R  14 
T2=- < 2 . 69 7529E-2 > *B I R 15 
TH=T 1 +T 2- 1.24121 70 

D 1  =  C 5 . 878 1 09E- 1 >  +  < 4 . 4 66 1 4E-3 >  *TH  +  C -3 . 998582E-4 >  *TH 12 
D2= C -5 . 3629 1 E-6>  *TH 13+  < 5 . 5383 1 3E-7 > *TH 14 
D3= C -7 . 405356E-9) *TH 15 
PD= ( D 1 +D2+D3  >  *  1 . 08080 

C1-C7. 065830E-4 >  *TH- C 3 . 83239E-4 >  *TH 12+ < 4 . 834588E-5 > *TH 13 
C2-C-1. 43039E-6  >  *TH  T4+  C  2 . 6 1 2 1 08E-8  >  *TH  15-  C  2 . 31 3786E-10>*TH 16 
CP=C 1 +C2 

MV* ABS C < P  < 5  > -PMOV  >  *2/ C RHT*PD > > : MV=SQR C MV > 

BT=P < 5  >  +CP*RHT*MV#MV/2 
I F  P <  1  >  >=P  C  3 >  ANDP f,'  4  >  >=P  <  2 >  THEN  QU=  1 
IF  PCI) >=P C 3 > ANDP C 2 > >=P C 4 >  THEN  QU=2 
IF  PCI) <=P < 3 > ANDP C 2 > >=P C 4 >  THEN  QU=3 
IF  PCD <=F C 3 > ANDP C 4 > >=P C 2 >  THEN  QU=4 
PF-BT/P1 

PRINTttl ,  "PI  " PC  1  > "PA"  ,  "QUADRANT" ; QIJ,  "  [THETA!  " ;  TH;  "DEGRES"  ■'  PRINTttl 

PRINTttl,  "P2";PC2>;  "PA",  "  IF  1 3  " ;  F I ;  "DEGRES"  ••  PRINTttl 

PR I NT# 1 , " P3 " ; P  C  3  > ; " PA " , " PD " ; PD , " VELOC I  TV " ; MV ; " M/S "  = PR I NT# 1 

PR  I  NT#  1 , "  P4 " ;  P  C  4  > "  PA " ,  "  CP " ;  CP , "  PR .  TOT .  LOC  " ;  BT ; "  PA "  =  PR  I  NT#  1 

PRINTttl ,  "P5" ; PCS) ;  "PA" ,  "PMOV" ; PMOV;  "PA" ,  "PP" ;  PF'  =  PRINTttl  : PRINTttl : PRINT# 

GO  TO  1850  , 

T 1 =2 1 . 53083#B I R- 1 3 . 94683*3 I R 12 

T2=8 . 953023*B I R 13- C  3 . 248728  >  *B I R#B I R*E I R*E I R 

T3=  C  5 . 903704E- 1 ) *B I R 15- C  4 . 1 4847E-2 >  *B I R 16 


.  -■w  r-  > 


55  -• 


1620  Dl-<5.  909329E- 1  >  +  68 . 1 934  rtJfc  -  4  J  *  I  H~  <.  6 .  rYVbloh-O)* i HTL 

1621  D2=<-5. 9604 1 5E-6) *TH 13- < 1 . 537824E-7)*THT4  H  ^ 

1622  B3=a.653353E-8>*TH15-62.34768E-10>*TH16  1 

1 623  PD= < D 1 +D2+D3 >  *  1 .  O 

1625  Cl=C-3. 25646E-3 ) *TH+  <  1 . 213657E-3>*TH12-<  1 . 436869E-4)*TH13 

1626  C2=  < 8 . 768080E-6 > *TH  14- < 2 . 1 72589E-7 ) *TH 15+  <  1 . 907035E-9>*7H  16 

1629  CF‘=C  1 +C2+3 . 777536E-4 

1630  GO  TO  1525 

1650  IF  DD<=0 . 0 1  THEM  1635 

1660  Tl=17. 266941EIR-8. 755624+EIRT2 

1 66 1  T2-5 . 74 8989 *B I R  T3-2 . 1 487921B I R  14 

1 662  T3=  < 3 . 965276E- 1 > *E I R 15- < 2 . 83 1 6G3E-2  > *B I R 16 

1663  TH=T 1 +T2+T3+0. 4941921 

1670  Dl=(5. 952E-1 )+<7. 559773E-5)*TH+6 1 . 564457E-4)*TH12 

1671  D2=*:~4.  3 4 5 3 9 5 E - 5  >  *  T  H 1 3 + (  2 .  1S2540E-6)*TH14 

1 672  B3=< -4 . 33095E-S  >  *TH 15+  C 3 . 24674 1 E- 1 0  > *TH 16 

1 673  PE=  <  D 1 +B2+D3 >  *1.0 

1 675  C 1  =  < -2 . S36S7E-3 ) #TH+  < 6 . 1 99820E-4 > #TH 12- < 3 . 45537SE-5 > #TH 13 

1676  02=  < 1 . 651 172E-6>*THt4-<2. 517231E-8)*TH 15+C7. 331682E-1 1 >*TH16 

1 677  CF-C 1 +02+ 1 . 402665E-3 
1680  GO  TO  1525 

1685  T 1  =  1 2 . 12255*£IR+C2.  32958 1 > * < B I R 12 > -3 . 229866*B I R 13+  < 1 . 2840€8)*BIRT4 

1686  T2=<-2. 268253E-1 )#BIR15+< 1 . 48998E-2)*£IR16 

1 687  TH=T 1 +T2-9 . 83739 1 4 

1 700  B 1  =  < 5 . S33236E- 1 , +  < 2 . 877764E-3 > #TH+<-l. 635622E-4 > *TH 12 

1701  B2=<-1 . 4225 1 6E-5  >  *TH  T3+ < 6 . 928993E-7)*TH14 

1 702  B3=  < -8 . 853089E-9 > *TH 15 

1 703  PD= <  D 1 +D2+D3 > * 1 . 80 

1705  01  =  (-8. 1 30340E-4 > *TH-  < 2 . @66893E-5>#TH12+<2. 829438E-5)*TH13 

1706  02= ( -4 . 975929E-7 ) *TH 14+ <  2 .  329837E-9>*TH 15 

1 707  CPC  1  +02+7 . 53634 1 E-4 
1710  GO  TO  1525 

1850  PRINT  "ACQUISITION  SUSPENDED .  FOR  NEXT  POINT,  TVF'E  1" 

I860  INPUT  KK 

1870  IF  KK=1  THEN  1310 

1900  SVS ( 3 1 362  > 

1918  PR  I  NT  "  fiCQIJ  I S I T I  ON  SUSPENDED .  TO  RESUME ,  TVF'E  1 " 


1920 

INPUT 

KK 

1 939 

IF  KK 

= 

1  THEN 

110 

I960 

GO  TO 

220G 

2200 

END 

5000 

X 1  =0  : 

V 

1=0  9 

XV 

=0 : X2=0 

5010 

FORJ= 

1T0NP 

5020 

X1=X1 

+ 

X<J> 

■  1 1 
■  T 

1=V1+VC 

5030 

XV=KV 

+ 

X<  J): 

*v 

CJ) 

5040 

V”'— VC* 

+ 

X  ( -J  > 

12 

:  V2=V2+' 

5050 

BN=X2 

- 

X1*X 

1/ 

NP 

5060 

A1CX 

1 1 
T 

-XI*' 

VI 

XNP ) XDN 

5070 

A0=<  V 

1 

-Al* 

XI 

)XNP 

5080 

R2=A1*A1*C 

-XI 1X1/1 

5890 

R2=SQ 

R 

(.  AES 

(R 

2)  .!> 

5100  RETURN 

8000  I  NF'UT  "  NUMBER  CHANNEL  " ;  N 
8020  POKE  40970,0 
8030  POKE  40971,0 

8040  H=PEEK< 40972 > : IF  NOTA  AND  123  THEN 
8050  A=PEEK C 40973) : IF  NOTA  AND  128  THEN 
8060  A=A-248 : GO  TO  8080 
8070  A=A+8 

8080  A=A*256+PEEK < 40974 ) -2048 
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FIG.  3D  -  CP  VERSUS  e  CALIBRATION  CURVES  FOR  FIVE  HOLE  PROBE 
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FIG.  9C  -  PD  VERSUS  e  CALIBRATION  CURVE  FOR  FIVE  HOLE  PROBE 
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FIG.  10  -  CALCULATED  ROLL  ANGLE  VERSUS  MEASURED  ROLL  ANGLE 
USING  FIVE  HOLE  PROBE 
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FIG.  M  -  FLOW  FIELD  EXISTING  AT  A  CLEAN  TEST  SECTION 
STATION  OF  THE  L-2A 
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THREE  VIEW  DRAWING  OF  HERCULES  C-130  AIRCRAFT  AND  FUSELAGE  STATIONS  WHERE 
PRESSURES  AND  FLOW  FIELD  MEASUREMENTS  WERE  TAKEN 
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FIG.  23  -  FLOW  FIELD  MAP  OF  FUSELAGE  STATION  2  ALONG 
BOATTAIL  OF  C-130 


*  *  i 


\  \  \  \  \  \  v 

UUUA 

\  \  \  \ 


\  \  " 

NX'' 


\  ' 


\  ^ 


.X  / 


// 
/  / 
/  / 


Sift 


t  t 


A  *  / 


A  •  * 


/  *  * 


\  *  A 


«  9 


/a  *  r  i  t- 


i  I 

/  1 

1  \ 

1  \ 

A  A 

A  \ 

*  r 

\  *t. 


FIG.  2A  -  FLOW  FIELD  HAP  OF  FUSELAGE  STATION  3  ALONG 
BOATTAIL  OF  C- 130 
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FIG.  25  -  FLOW  FIELD  MAP  OF  FUSELAGE  STATION  4  ALONG 
BOATTAIL  OF  C-130 
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FIG.  26A  -  FLOW  FIELD  MAP  OF  FUSELAGE  STATION  5  ALONG 
BOATTAIL  OF  C- 130 
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FIB.  26B  -  ISOLINES  OF  TOTAL  PRESSURE/DYNAMIC  PRESSURE 
ON  FUSELAGE  STATION  5 
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FIG.  27B  -  ISOLINES  OF  TOTAL  PRESSURE/DYNAMIC  PRESSURE  ON  FUS.  STA.  6 

-  CHARACTERISTIC  ISOLD'S  OF  PRESSURE  LOSS  DUE  TO  WAKE  INDUCED  BY 
AIRCRAFT  BOATTAIL. 

-  MEASUREMENTS  TAKEN  AT  1  FUSELAGE  DIAMETER  DOWNSTREAM  OF  FUSELAGE 
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FIG.  27C  -  ISOLINES  OF  CROSS  FLOW  KINETIC  ENERGY  AT 
FUSELAGE  STATION  6 
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FIG.  35  -  STRAKES  CONFIGURATIONS  (H  =  0,  0.5,  1.0) 
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FIG,  38B  -  ISOLINES  OF  DIMENSIONLESS  CROSS  FLOW  KINETIC  ENERGY 
(DCFKE)  DOWNSTREAM  OF  H  =  0,  ap  =  -10  THIN  STRAKE 
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FIG,  m  -  ISOLINES  OF  P„  DOWNSTREAM  OF  H  =  0,  «„  =  15  THIN  STRAKE 
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FIG.  VC  -  ISOLINES  OF  DCFKE  DOWNSTREAM  OF  H  =  1.0,  %  =  15  THIN  STRAKE 
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FIG.  51C  -  ISOLINES  OF  DCFKE  DOWNSTREAM  OF  H  =  0.5,  c*R  =  15  THICK  STRAKE 
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FIG.  S5B  -  PRESSURE  DISTRIBUTION 


FIG.  56A  -  SCHEMATIC  OF  BOATTAIL  INDUCED  WAKE  PATH 
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